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Abstract 


Laminated composite structures exhibit a number of different 
failure modes. These include fiber-matrix debonding within individual 
layers, delamination or separation of the layers, transverse cracks 
through one or more layers and fiber fracture. These failures are 
influenced by enviromental conditions. Thermal and moisture 
conditions are significant factors in interlaminar delamination as 
well as the other modes of failures. 

A simple delamination analysis method is presented here. It is 
based on a shear -type deformation theory and includes hygro thermal 
effects. These enviromental conditions are applied to the strain 
energy release rate and interlaminar shear stresses. 

The method is applied to mixed mode edge delamination specimens 
made of T300/5208 graphite/epoxy material. Residual thermal and 
moisture stresses significantly influenced the strain energy release 
rate and interlaminar stresses. Both experienced large increases when 
thermal conditions were added to the mechanical strains. These 
effects were alleviated when moisture stresses were included. Thermal 
effects on the interlaminar shear stress and total energy release rate 
were totally alleviated for the same specific moisture content. 
Moreover, this value of moisture content was not significantly 
affected by the stacking sequence for the laminates considered. 
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Introduction 


Composites have been used in the aircraft industry since 1969. 
One aspect of concern for using composites in structures is separation 
of plies or delamination. This occurs in regions of stress raisers 
such as holes, ply terminations, cut-outs and free edges. 
Delamination along the free edge of laminates have been observed 
during testing and service. The presence of delamination, initiated 
by interlaminar stresses, causes redistribution of the stresses among 
the plies in a laminate. Thus, it usually results in a reduction of 
stiffness and strength. 

Figure 1 shows delamination in a rotor hub made of S2/SP250 

glass/epoxy. The specimen has delaminated in two places that can be 

depicted by the dark lines. Figure 2 shows delamination in a single 

cracked-lap-shear test specimen made of AS4/3502 graphite/epoxy [ 1 } . 

The specimen layup is [+45/0/90] quasi-isotropic with 8 plies in the 

6S 

strap and 40 plies in the lap. The tests were performed on a 
displacement controlled machine. Fiber glass tabs were attached to 
the specimen ends. Multiple, isolated free edge delaminations occur 
in the neighborhood of the lap/strap junction and the tab. Figure 3 
illustrates an I-beam section made of C3000/5225 graphite/epoxy woven 
cloth in the post-buckled regime. Free edge delamination depicted in 
the flanges precipitated the final failure in the specimen. Figure 4 
shows how delamination can take place in single cracked- lap - shear 
specimens subjected to compressive and tensile loading. Specimens A 
and B delaminated under compressive loading while C experienced a 
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tensile loading. These examples illustrate the importance of 
investigating delamination problems in composites. 

Thermal residual and moisture effects on a composite are 
practical enviromental conditions. Determining the response of these 
conditions on interlaminar stresses and energy release rates in 
laminated composites is the primary objective of this work. 

Delamination analysis can be based on two approaches. They are 
the strain energy release rate and interlaminar stresses. The 
interlaminar stresses are due to Poisson's ratio mismatch and 
difference in the coefficients of thermal and moisture expansion 
between plies. Delamination occurs when these stresses reach the 
interlaminar strength of the material. An alternative approach is 
based on the actual process of fracture rather than the strength 
concept. Delamination can propagate when the strain energy release 
rate at the crack front is sufficient to overcome the material's 
fracture resistance or toughness. 

The strain energy release rate can be obtained for three 
particular modes of crack action. These three modes are known as Mode 
I or opening mode, Mode II or forward shearing mode and Mode III or 
tearing mode. Several failure laws are formulated in terms of these 
modes { 2 ] . 

A simple analysis predicting interlaminar shear stresses and total 
energy release rate with the influence of thermal and moisture effects 
is developed. This simple approach is useful in understanding the 
basic mechanics of the problem and predicting the factors controlling 
the behavior. The method is directed to the needs of a practical 
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design environment. It is not intended to compete with large-scale 
numerical approaches, but rather to serve as the means for selecting 
and screening candidate configurations and providing trend 
information. Simple codes for a desktop computer have been created to 
analyze laminate configurations. 
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Literature Summary 

A historical discussion of previously developed work for 
predicting interlaminar stresses and energy release rates is presented 
to establish a basis for the proposed model and to permit the present 
work to be placed in proper perspective. 

Earlier analyses have reflected the prediction of interlaminar 
stress and energy release rate without hygrothermal conditions. 
0'Brien[3,4] investigated delamination onset and growth in 
graphite/epoxy laminates under uniform extension. A simple expression 
was developed for the total energy release rate based on classical 
lamination theory. Whitney and Knight [5] used classical laminated 
plate theory to develop an edge delamination specimen analysis. This 
work was limited to Mode I behavior. 

An analysis based on a shear deformation theory and a sublaminate 
formulation [6] was developed by Armanios and Rehfieldf 7 , 8 ] . This 
method provides good estimates for the interlaminar shear stresses. 
Energy release rate components are estimated based on these stresses. 
However, this method does not provide reliable estimates of peel 
stress since thickness strain is neglected. This analysis was limited 
to mechanical strain only. 

0 / Brien[9] modified his analysis to include thermal and moisture 
conditions. The influence of thermal effects was considered by 
Whitney ( 10] . 

The work of Reference 9 was based on a classical laminated plate 
theory. It was applied to mixed-mode edge delamination specimens. 
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The results were limited to strain energy release rate. Finite 
element modeling was used to determine the strain energy release rate 


energy release rate due to thermal effects. It decreased with the 

addition of moisture considerations . For a T300/5208 graphite/epoxy 

laminate with [±30/+3 0/9 0/93] layup, the thermal effect increased 

I s 

the total energy release rate by 170 percent when compared to 
mechanical loading alone. However, a specific moisture level of 0.75 
percent completely alleviated this increase. In calculating the total 
strain energy he showed that bending and coupling effects became 
important at high levels of moisture content. 

In Reference 10 a higher-order plate theory with transverse 
normal strain effects was developed. Peel as well as interlaminar 
shear stresses could be predicted by this method. The thermal 

influence on total energy release rate and interlaminar stresses was 
investigated using a Mode I specimen. Residual thermal effects showed 
a significant influence on the stresses and release rates. For a 
graphite/epoxy laminate of [O 3 , 9 O 3 ] layup , thermal effects increased 
the maximum peel stress by a factor of 2.7 over that of pure 


mechanical strains. 


In the present work both thermal and moisture influences are 
studied in a mixed-mode delamination specimen. The analysis includes 
total energy release rate as well as interlaminar stresses. 
Similarities between the interlaminar stesses and total energy release 
predictions with hygrothermal effects is investigated. 

In the subsequent sections, the analytical approach is developed. 
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The method is then applied to six graphite/epoxy laminates. A 
discussion of the hygrothermal effects on interlaminar shear stress 
and total energy release rate predictions is provided. 
Recommendations for further investigations are proposed. Appendices 
are included for completeness . The first provides detailed 
expressions of the governing equations. Appendix II defines the 
hygrothermal expressions and their use in the analysis. The last 
appendix shows a listing of the program used and sample output. 
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Analytical Approach 


Overview 

The sublaminate modeling approach describes the essential 
features of the laminate behavior in a simple way. A free edge 
delamination specimen is shown in Figure 5. A uniform strain , €, is 
applied in the axial direction. From symmetry, only one quarter of the 
specimen is considered. In Figure 6, the specimen is modeled as if it 
were composed of four distinct sublaminates. Sublaminates 2 and 3 
represent the group of plies above and below the crack, respectively 
in the cracked portion of the laminate, while sublarainates 1 and 0 
denote the same group of plies in the uncracked portion of the 
laminate . 

The use of sublaminates -- groups of plies that are conveniently 
treated as laminated units -- simplifies the analysis considerably. 
This approach is applied with confidence when the characteristic 
length of the response is large compared to the individual sublaminate 
thickness [6] . This sublaminate modeling approach has been verified in 
Reference 7 by comparison with a ply-by-ply finite element solution. 
These sublaminates are connected by enforcing the proper continuity 
conditions on stresses and displacements at their interfaces. 

Displacement fields within each sublaminate are defined as: 

u - x€ + U(y) + z(3 x (y) 

v - V(y) + z/? y (y) (1) 


w - W(y) 
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where u,v, and w denote the displacements relative to the x, y, and z 
axes, respectively. Shear deformation is recognized through the 
rotations JS X and £y. The governing equations for each sublaminate are 
derived using a virtual work approach. The derivation of the 
governing equations used in the development appears in Appendix I. The 
derivation is an extension of the work of Reference 8 with 

hygro thermal effects included. 

The constituitive relationships in terras of these force and moment 
resultants can be written as 

Nj-A^j +B ik »c k -N^ M (i,j,k - 1,2,6) 

Mj - Bjj € j + D jk < k - Mf (i,j,k - 1,2,6) (2) 

Qi - Ajj (i,j - 4,5) 

where the subscripts x, y, z, yz, xz, and xy are replaced by the 
subscripts 1-6 respectively. The force and moment resultants are 
denoted by N t and M ( - , respectively. Non-mechanical forces and moments 
resulting from the hygrothermal effects are labeled with a superscript 
NM. They are defined as: 

NM NM 

(Nj , M ) -/ Qij(ltZ) (Of. AT + b. C) dz (3) 

-h/2 J J 

The swelling coefficient is denoted by b j in Equation (3), the 
thermal coefficient by Oi- . The change in temperature is denoted bv 
AT and moisture weight gain by C. 
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The elastic stiffnesses A jj , By , and Dy are defined in terms of 
the sublaminate reduced stiffness Qy for a plane stress situation. 
These bear the classical definition. 

h/2 

< A ij • B ij • D ij > -jT /2 (l.z.z 2 ) Q j j dz (4) 

The equilibrium equations are: 

^xy,y + " t lx^ “ 0 

Ny,y + (t 2y - t ly ) - 0 


Qy,y + (P 2 

V 

- 0 



(5) 

**xy, y * Q x + 

h/2 

< c 2x 

+ t lx ) 

- 0 


M y,y ■ Qy + 

h/2 


+ t^ x ) 

- 0 


where t 2x , t 2 y, p 2 and 

c lx 1 

t iy* 

p^ denote the 

interlaminar s tress 


components in the x-z, y-z and z directions at the sublaminate upper 
and lower surfaces, respectively. These stress components appear in 
Figure 7. 

The displacements, resultant forces and moments, and interlaminar 
shear stresses in each sublarainate is governed by the displacement 
distribution (1), constituitive (2), and equilibrium (5) equations. 
These equations are applied to each sublaminate. The variables 


associated 

with 

each 

sublaminate are coupled through 

the 

continuity 

requirements 

at 

their 

interfaces . 

Enforcement of 

the 

boundary 

conditions 

lead 

to a 

solution for 

these variables. 

This 

procedure 
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discussed in general terms above is applied to the analysis of the 
edge delamination specimen shown in Figure 2 in the following sections. 

The response associated with sublaminates 1 and 0 shown in Figure 
2 is coupled through the continuity conditions at their common 
interface. The situation is different with sublaminate 2 and 3 where 
the continuity conditions are relaxed due to the presence of the 
crack. 

Uncracked Region: Sublaminates 0 and 1 

From symmetry conditions at the sublaminate bottom surface the 
shear stresses are zero. Interlaminar stresses at the top surface of 
sublaminate 0 are equal to those on the bottom of sublarainate 1 . 
Substituting these conditions into the equilibrium and constitutive 
relations and enforcing continuity of displacements at their common 
interface yields a homogeneous system of ordinary differential 

equations. These can be expressed in terms of the sublaminate 
rotations /J x and )3y. Assuming an exponential solution of the form 

(0ix*^ox »0iy >0oy ) - (0ix,£ox »0iy »^oy ) (6) 

results in a characteristic equation of the form 

EgS 8 + E 6 s 6 + E 4 s 4 + E 2 s 2 + E q - 0 (7) 

Parameter Eq depends only on the stiffness coefficients A^ 4 , A^ 
and A 45 for both sublaminates while Eg is predominantly influenced by 
the bending and coupling coefficients D - j and B,j . 


Thus, its 
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numerical value can be orders of magnitude smaller than the other 
coefficients. This results in the presence of a boundary zone in the 
response . 

The characteristic roots controlling the behavior of the laminate 
are determined from Equation 7 which has a closed- form solution. 

Crack Region of the Laminate: Sublaminates 2 and 3 

With this group of laminates, there are free surfaces at both the 
top and bottom of sublaminates 0 and 1 respectively. This is due to 
the presence of the crack. With the crack dividing the sub lamina te s , 
continuity conditions are not enforced at the boundary interface. 
This results in zero shear stresses at the surfaces of each 
sublaminate. Thus, the equilibrium and constitutive relations combine 
to produce a second order differential equation in terms of the 
sublaminate rotations #2X f° r sublaminate 2 and /? 3 y for sublaminate 3. 

'Interlaminar Stresses 

The arbitrary constants that are obtained from the eighth degree 
polynomial are determined by enforcing the stress free boundary 
conditions at the free edges of sublaminates 2 and 3, and the 
continuity of force, moment, displacement and rotations between 
sublaminates 0 and 3, as well as between 1 and 2. This yields the 
following expression for the interlaminar shear stresses. 

C X - "XYJ.Y- °j "' S,y (8) 

C Y - %,Y - T J 


( j “ 1 - 4 ) 


(9) 
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Parameters Tj and Gj represent the amplitude of the response. 

Energy Release Rate 

A complete formulation of the strain energy release rate appears 
in Appendix II. The total energy release rate can be determined by 
considering the work done by external forces. 

G 1/2 * dW/da (10) 

The total energy release rate is denoted by G and the crack 
length by a. The work done is defined as 

w - l/2 -C Nxi^mi d y < U > 

where subscript i denotes the sublaminate being referenced. The term 
represents the mechanical strain in each ply. This is defined as the 
difference between the total strain and the strain corresponding to 
free expansion for each ply. This strain is estimated by using a 
procedure similar to Reference 5. However, in Reference 5 the free 
expansion strain was determined by considering groups of plies in the 
cracked and uncracked regions of a Mode I edge delamination specimen. 
This approach is valid for a limited class of laminates. A 
ply-by-ply analysis rather than a sublaminate modeling should 
be used. In the following analysis, free expansion strains are 
determined on a ply-by-ply basis. 

From the symmetric conditions that exist in the uncracked section 
of the laminate, there exist no curvature. In the cracked portion. 
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the moment about the y-axis is assumed to be zero. Using both of 
these boundary conditions in Equation (2) yields the following. 

NM 

N Xl“ A ll € + A 12 € y + A 16 € z * N xl ^ ^ 

»,NM 

N xl“ A 11 €+ a 12 € Y + a 16 €z + B 12 K y" n x2 


Strain components € and G in Equation 12 are expressed in terms 

y z 

of the applied strain by 


„ ,NM 

Gy- C v € + C v 


(13) 


€ Z- C u €+ c 


NM 

u 


NM NM ^ . . 

The terms . C v , C u , C v and C u are functions of the extensional 

stiffness components of sublaminates 1 and 0. 


Using these expressions, Equation 12 can be re-written in the 

form. 


._k . k .^k _k . 

N xc - ( E c « c + T c ) 


(14) 


- ( E k € k + T k ) 
^XU ^ u X U * 


k k k k 

Parameters E c , T c , E u and T u are defined in Appendix II. 
Superscript k denotes the individual ply. Subscripts u and c 
represent the uncracked and cracked portions, respectively. The 
non-mechanical strain in each ply corresponding to a state of free 



17 


expansion is obtained by allowing the stress in each ply to vanish. 
This yields the following. 


k 

«u 


_k . 

- T u / 


E k u 


(15) 


k 


-Tc / Ec 


The strain corresponding to free expansion of the entire laminate 
is obtained by letting the resultant force vanish. The non-mechanical 
strain is 


€ NM - (t* - (tJ - T* ) 2a/b ]/| eS - ( eJ - E* ) 2a/b} (16) 


4c sic 4^ a k if k 

The terms Tj\ T c , E u and E* represent the summation of T u , T Q , E u 

k 

and E c over their respective sublaminates . These strain definitions 
for the effects of moisture and temperature can now be used in the 
general expression for the strain energy. The strain field is altered 
to represent the hygrothermal effects. The total strain for a 
sublaminate is expressed as: 


€ + € 


NM 


(17) 


The strain energy expression is given below showing the use of 
Equations (13) and (17). 

w - I ( € T + t£) (<= T - ) + ( Eu e T + t£ ) ( € T - ) j 


(18) 



Substituting this 
energy release rate 
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into Equation (11) yields the total strain 
per 


crack. 
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Results and Discussion 


Benchmark Study 


The analytical model is applied to the edge delamination specimen 
shown in Figure 5. The material considered is T300/5208 
graphite/epoxy . Its mechanical properties are listed in Table I. The 
coefficients of swelling and thermal expansions are also stated. The 
geometry of the specimen is given in Table II. Cure temperature for 
this material is 350° F. The ambient operating temperature is 70° F. 
The moisture level for all cases was allowed to vary from 0 to 1.2 
percent of the laminate weight. This moisture level reflects feasible 
conditions. The mechanical strain is taken as 0.00254. This 
particular value of strain was taken from test experiments [ 9 ] . It is 
considered a practical value for the material. 

Six laminates have been analyzed. They can be divided into two 
groups. The first group is composed of laminates [35/ -35/0^90^, 

I 35 /o / - 35 f 01 >- and t° / 35 / - 35 't 90 V la ” lna “ s are pr ° n<i co 

delaminate at the interfaces indicated by the arrow[9] . The Mode III 
in these laminates is negligible. This is due to the fact that 
relative sliding at the crack front and the interlaminar shear stress 
in the x-z direction, T Xz , is neglegible. The second group of 

laminates is [ 30/-60'75/15 L [ -35/55/10/-80]. and [ 35/80/- 55/- 10 ] . In 

f S' | s I s 

these laminates Mode I, Mode II, and Mode III are finite. The Mode 


III strain energy release rate component due to mechanical loading in 
these laminates are significantly large, ranging from 60 to 90 
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percent^ 7] . 


Interlaminar Shear Stresses 


The interlaminar stress at the delaraination interface appear 
in Figure 8-10, for the first group of laminates. The interlaminar 
shear stresses T XZ and Ty z for the second group of laminates appear in 
Figures 11 - 13. The labels M, M+T, and M+T+H stand for mechanical, 


mechanical and thermal, and mechanical, thermal and moisture 
respectively. 

The boundary layer of decay for all laminates ranged from 0.85 to 
0.93 of the laminate semi-width. In this context the boundary layer 
decay length is defined as the distance where the stress decays to 
5 percent of its maximum value. -The stress boundary zone is not 
significantly influenced by the environmental conditions. 

The magnitude of shear stress however showed a strong dependency 
on thermal and moisture conditions. At the delamination front, the 
ratio of stress with thermal effects as compared to pure mechanical 
loading ranged from 3.22 to 3.36 for the first group of laminates. 
This maximum was experienced at the crack tip. For the laminates 
where Mode III was present, this ratio ranged from 4.16 to 5.23 for 
Ty Z . The shear in the x-z direction showed a ratio of 1.4 to 2.16 for 
the maximum stresses. The maximum Ty Z stress for the second group of 
laminates was experienced at the crack front. However, the maximum 
T'xZ stress occured slightly ahead of the crack. This can be seen 
In Figures 11-13. 



21 


The addition of moisture alleviated the thermal effect. A 
moisture content of 0.4 has reduced the stress of thermal influence by 
approximatly 40 percent as compared to thermal influences alone. This 
trend is similar to the results of Reference 9. 

Numerical values of interlaminar shear stress at various moisture 
levels are provided in the sample output of Appendix III. 
Interlaminar shear stresses show numerical decrease with increase of 
moisture levels. 

Energy Release Rate 

The hygrothermal effect on total energy release rate appears in 
Figures 14 -15. The hygrothermal effects on total energy release rate 
show a similar trend to that of interlaminar shear stresses. Residual 
thermal stress tends to increase- total energy release rate while 
residual moisture alleviates this effect. The figures show that for 
total alleviation of the thermal effect, the specific moisture content 
ranges between 0.70 and 0.77 for all laminates. This indicates there 
exist a weak dependency on the stacking sequence. 

The effects of thermal conditions alone on the energy release 
rate does not correspond to the same numerical value as the 
interlaminar shear stresses. The total energy release rate of layups 
where Mode III was negligible showed a ratio of 5.1 for mechanical and 
thermal compared to mechanical conditions only. For the laminates 
where Mode III is finite, this ratio varied from 1.6 in the 

( 30/ - 60/75/- 15 ) layup to 3.37 in the [ 35/80/- 55/- 10 1 layup . 

T s I s 
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The total energy release rate in the first group of laminates is 
approximately the same for mechanical loading as shown in Figure 14. 
The influence of thermal and moisture does not appear to alter this 
trend. The energy release rate for the [ 35/- 35/j0/90] s laminate is 
indistinguishable from the [0/35^-35/90^ layup. The rate of 
alleviation due to moisture is the same for the three laminates. This 
is in contrast with the alleviation rate of the laminates where Mode 
III is finite as shown in Figure 15. For this class of laminates, the 
rate of alleviation due to moisture is different for each laminate. 

In some of the laminates, the rate of alleviation is not 
constant. There is a steep gradient in the rate of alleviation until 
the moisture content approaches the totally alleviated state. After 
such moisture content, the decrease in total energy release rate with 
respect to moisture addition is not as significant. 

■ It is worth noting there is a similarity between the strain 
energy release rate prediction and the interlaminar stresses for the 
totally alleviated state. This is shown in Figure 16 for a 
[ -35/55^10/-80] s layup. The specific moisture percent producing 
complete alleviation of the total energy release rate from the thermal 
effect is 0.76 as seen in Figure 15. The interlaminar shear stress 
distribution corresponding to this level of moisture is 
indistinguishable from the mechanical loading alone. The same 
conclusion was reached studying the other laminates. 



Conclusions 


A simple analysis was developed that predicted the influence of 
thermal and moisture effects on the interlaminar shear stresses and 
strain energy release rate. The analysis was applied to six 

mixed-mode edge delamination specimens. The results provide several 
significant findings. 

1. Residual thermal strain has a significant influence on the 
interlaminar shear stress and total energy release rate. 

The interlaminar stress and total energy release rate 
increased by 330 and 510 percent respectively over that of 
pure mechanical loading. 

2. Moisture tends to alleviate the thermal effect for both the 
interlaminar stress and energy release rate. At a specific 
moisture content of approximately 0.75 percent, the thermal 
influence is totally alleviated. 

3. The moisture content for total alleviation found from the 

total energy release rate analysis also produced an 

* 

interlaminar stress distribution similar to pure mechanical 
loading conditions. 

The first two findings are in agreement with the results of 
previous investigators. The third finding is new. It establishes a 
similarity in behavior between a delamination analysis expressed in 
terms of the energy release rate and the strength approach expressed 
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by the interlaminar stresses. 

These findings point to new directions for further inquiry. 
These are discussed in the following section. 
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Recommendations 


The thermal effects on the laminates showed a large increase in 
both the interlaminar shear stresses and strain energy release rate. 
The analysis should be supplemented with experimental tests to verify 
the result. Several fracture laws are expressed in terms of the 
strain energy release components, as well as the total strain energy. 
Further analysis should include predictions of these components in the 
presence of hygrothermal conditions. 

Throughout this work the temperature is assumed to be uniform 
through the thickness of the laminate. The same is true with the 
moisture. An approach corresponding to a practical environment method 
should account for temperature and moisture gradients in the laminate. 
In this situation, the hygrothermal gradients through the thickness 
may create an unbalance effect in an originally balanced construction. 
This consideration is of significant importance in aerospace 
structural components subjected to a large temperature difference 
between the upper and lower surface. 

The loading considered here is uniaxial extension. However, it 
is known that the load transfer points are not always in the plane of 
the laminate. Therefore, investigating laminate response under 
combined loads is' of great practical importance. It is recommended 
that bending, torsion as well as their combined effect be addressed. 

Findings by previous investigators suggested that delamination 
behavior in laminates subjected to fatigue loading follows static 
loading conditions. Further work is needed to investigate the 
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influence of hygrothermal conditions on the delamination of laminates 
under fatigue loading. 

Finally, the present analysis is applied to the mixed-mode edge 
delamination specimen. Extension of this work to other specimens such 
as the single- and double -crack- lap shear and the Mode II edge notch 
flexure specimen is recommended. 

When accomplished, these recommendations, together with the 
present research will provide a better understanding of the 
delamination problem in composites. Consequently, this will enable 
predicting, managing and ultimately preventing interlaminar fracture 
in laminated composites. 
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TABLE I - T300/ 5208 GRAPHITE/ EPOXY PROPERTIES 

£ 12 = 18.7 MS I 

E 22 = 1 .23 MSI 

G 12“ °* 8 3 2 MSI 

Poisson Ratio *= 0.292 

Swelling Coefficients of the Material direction: 
b ( 1 -di rection) ■ 0 
b ( 2 -direction) - 556 O b e/ %w eight 

Thermal Coefficients of the Material direction: 
tt(l-di rection) ■ -23 Ue/*f 
Qt(2-di recti on) • 1b«9uc/*F 


TABLE II - GEOMETRIC DIMENSIONS OF SPECIMEN 


Ply thickness « 0.005b inch 
Width = 1.5 inch 

Crack length ■ 6 x ply thickness = 0.032b inch 
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FIGURES 1-16 
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FIGURE 1 - GLASS /EPOXY ROTOR HUB AFTER DELAMINATION 
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r~ «3E IS 

QUALITY 





FIGURE 2 - GRAPHITE/ EPOXY SINGLE CRACK- LAP- SHEAR SPECIMEN 
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GRAPHfTH/ EPOXY POST- 
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FIGURE 5 - FREE EDGE DELAMINATION SPECIMEN 
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FIGURE 8 - COMPARISON OF INTERLAMINAR STRESS 




FIGURE 0 - COMPARISON OF INTERLAMINAR STRESS 
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FIGURE 10 - COMPARISON OF INTERLAMINAR STRESS 




FIGURE 11 


COMPARISON OF INTERLAMINAR STRESSES 
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FIGURE 12 
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FIGURE 13 - COMPARISON OF INTERLAMINAR STRESSE 



CNi-NI) / CN1-HD 


44 



FIGURE 14 - ENERGY RELEASE RATE DISTRIBUTION 
FOR LAMINATES WrTHOUT MODE ni 


i 
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FIGURE 16 - TOTAL ALLEVIATED STATE STRESS DISTRIBUTION 


APPENDIX I 
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Derivation of the Governing Equations 

In this Appendix the governing equations for the sublaminate shown in 
Figure 7 are derived using the principle of virtual work. 

Consider a sublarainate of thickness h. The origin of a cartesian 
coordinate system Is located within the central plane (x-y) with the z-axis 
being normal to this plane. The material of each ply Is assumed to possess a 
plane of elastic symmetry parallel to xy as shown in Figure 6 . 

Stress and moment resultants are given below. 


(N X . Ny, N X y , 

,h/2 

Qx» Qy) " J^x. CT y» r xy, r x 2 » r yz) 

dz 


(H X # My • M X y) 

rV 2 

” 1 (^x* °y> f xy) z ^ 

■*h/2 


(i-D 

Because of 

the existence of a plane of 

elastic 

symmetry, the 


constitutive relations are given by 


*x 


Cll 


c x 

°y 


c 12 C 2 2 SYM 


c y 

°z 


c 13 c 23 c 33 


c z 

r xy 


c 16 c 26 C 36 C 66 


Yxy 


r rz 

. 

C 44 

SYM 


lyz 

r xz 


C 45 

c 55 


Yxz 


where C^j are components of the anisotropic stiffness matrix and 7 X y, 7y Z and 
7 xz are engineering shear strains. 

The displacements are assumed to be of the form 


u - U(x,y) + z0 x (x,y) 
v - V(y) + z0 y (x # Y) 
w - W(x,y) 


d-3) 
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where u t v and w are the displacement components in the x t y and z directions, 
respectively . Equation (1-3) in conjunction with the strain-displacement 
relations of classical theory of elasticity leads to the following kinematic 
relations 


c xx “ u *x* Z ^x,x 
€yy — V,y + z Py 9 y 

c zz “ 0 

7xy U*y + z (6 x ,y+ Py 9 x) 

7xz “ £x ♦ W,x 

7yz ~ 0y + w .y < X “ 4 ) 


Substitute Equation (1-4) into Equation (1-2) and put the results into 
Equation (1-1). This yields the following constitutive relations: 




A 11 

a 12 

a 16 

B n 

b 12 

N y 


a 12 

a 22 

a 26 

b 12 

b 22 

^xy 

— 

a 16 

a 26 

a 66 

b 16 

b 26 

«x 


B 11 

b 12 

b 16 

Dll 

d 12 

H y 


b 12 

b 12 

b 26 

Dll 

d 12 

M xy 


b 16 

b 26 

b 66 

Dl6 

d 26 

1 

r 


r 




Qy 


A 44 

a 45 

Qx 


a 45 

a 55 


where 

f h/2 

(Aij , By, Djj) - J C tJ (1, z, z 


b 16 


D.x 


N x 

b 26 


V.y 


Ny 

b 66 


U,r+V, x 

- 

N X y 

Dl6 


^x,x 


«x 

D26 


^y.y 


My 

d 66 


0x, + 0r,x 


«xy 


fiy + W t y 

0X + 


(1-5) 


and the non-mechanical terms are defined in Appendix II. 

The variation of the strain energy due to virtual displacements 6 u, 6v 


i 


and 6 w is 
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i 


r* 


i 


SV 


- [ (°X 


Sc x + CyScy + a z Sc z + r X y ^7 X y + Ty Z Sy yz + r xz 6y xz )&V 


d-6) 


where ic x , Scy , Sc z . 6 7xy» $7xz are t ^ ie strains associated with the virtual 
displacements. Using Equations (1-3) and (1-1) then integrating through the 
thickness gives 


SV 


- | [N x x + Ny5V,y + N X y 5U,y + Q x S0 X + Qy ( 60y + SW.y) + H X S 0 X t x 


+My 60 


y.y + M xy 


£ £x f y 1 dA 


d-7) 


The variation of the work done by the external forces and by the 
surface tractions is 


SW •» I (n x 5U + ny SV + q6V + m x 80 x + niy 80 y) dA 
J A 

+ f («; «U n + N ns 6V S + M n 60 n + M ns *0 s )ds (1-8) 

J S 

where a bar denotes values on the boundary. Variables n and s are coordinates 
normal and tangential to the edge , and 

n x ” c 2x * c lx 


n y - t 2 y - t ly 

q - P2 - Pi 

m x " | ( c 2x + c lx) 

m y “ “ (t2y + tly) (1-9) 


where n x and ny can be regarded as effective distributed axial forces. Terms 
m x and my are effective distributed moments and q is an effective lateral 
pressure. 


From the principle of virtual work the equations of equilibrium and 
boundary conditions are determined from the Euler equations and boundary 
conditions of the variational equation. 

SV - 6W 




(I-lCn 
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Substitution of Equations (1-7) and (1-8) into Equation (1-10) leads the 
following equations of equilibrium: 

N x,x + N xy,y + n x - 0 

**xy,x + Ny,y + Hy - 0 
Qx,x + Qy,y + Q “ 0 
**x,x + **xy,y - Q x + m x - 0 

^xy , x ^y , y “ Qy + ®y ”0 (1-11) 

and one member of the following five products must be prescribed on the 
sublaminate edges 


^n ^n » N ns u s * ^n^n » ^ns @s an< ^ Qn ^ 


( 1 * 12 ) 


For the ED specimen under uniform extension, U(x,y) in Equation (1-3) 
is given by 


U(x,y) - U*(y) + x€ 


(1-13) 


and the response is a function of y and z coordinates only. For this case the 
equilibrium equations (I -11) take the form 


N xy,y + n x - 0 

Ny ^ y ^ Hy ™ 0 

Qy,y + q - 0 
^xy,y ‘ Qx + m x “ 0 

My,y - Qy + my — 0 (1-14) 


Substitution of the constitutive relations in Equation (1-5) into Equation (I- 
14) yields the following equilibrium equations in terras of kinematic 
variables . 



52 


A 66*-yy 


U* 


n x 

A 26*-yy A 22^yy SYM 


V 


n y 

0 0 -A^Lyy 


V 

- - 

-q 

^66^yy ®26^yy ~ A 45*-yy (^66^yy * A 55) 


fix 


m x 

®26^yy ®22^yy “^ 44^-7 ^26^-yy " A 45 (^22^yy ~ A 44) 

t 

A 

1 

niy 


( 1 - 15 ) 

where the operators 


Lyy - d2/dy 2 

Ly - d/dy 

From these governing equations the basis of the work presented in this 
paper is formed. Appendix II gives a detailed formulation of the 
hygr c thermal terms and the formulation of the total strain energy release 
rate and interlaminar stresses. 
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APPENDIX II 
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Hvgro thermal Effects on E dee Delamlnation 


The displacement field and constitutive relations governing the free 
edge ply separation were presented in Appendix I. The hygro thermal 
expressions, represented with the superscript NM for non-mechanical, are 
defined as follows 


< Hf ) 



(II-l) 


where 

aj — Coefficient of thermal expansion 
— Swelling coefficient 
T — Local temperature 
T r — Reference temperature 
C - Specific moisture concentration 
Qij - Reduced stiffness coefficient 


The terms aj and are transformed as second order tensors with the 
assumption of no thermal or swelling shear strain. 

The concept of eublaminates is used when enforcing the boundary 
conditions. 


Cracked Sublaminates 
Sublaminate 2: 

The boundary conditions for this sublaminate are expressed as: 

Ny2 - N xy2 - My2 - Qy2 - 0 (II-2) 

M xy2,y * Qx2 “ 0 

Using the first three conditions in the governing equations, one can 
express V 2 , U 2 and /?2y * n terms of ^2x to obtain: 



55 


42 


r ^ 

€ 

+ 

A 1 A 1 R 1 ‘ 

22 A 26 22 


v 2,y 


Nyl 



02x,y 


111 





A 16 B 66 


A^ a^ 

*26 A 66 A 26 


v 2,y 

- 

N xyl 

1 1 


1 1 1 





B 12 D26 


B 22 B 26 °22 


02y,y 


Mri 

W 4 


Sublaminate 3: 

The boundary conditions for this sublamlnate are given as: 
N y3 - N xy3 - 0 

M y3,y * Qy3 ** ® 

M xy3,y - Qx3 " 0 


(H-3) 


(II-4) 


These are used in a similar manner (as in sublaminate 2) to 




v 3,Y 

+ 

*26 


«3.Y 



A° R° R° 
x2 B 22 B 26 


.o _o _o 
A 16 B 26 B 66 


03, y 
03x,y 


N 


yo 


N 


xyo 


obtain 


- 0 


(II-5) 


These equations are then substituted back into the governing equations 
to obtain expressions for the force and moment resultants. They can be 
expressed in terms of the strain plus non-mechanical effects. 


UNCRACKED SUBUVHINATES 


Sublaminates 0 and 1: 


The boundary conditions of continuity at the interfaces must be 
satisfied. 


Ny Pb N xyl (0) - N y< £0>- N xyo <0)= 0 


(H-6) 


M X yl (0) - Mxy2(0) 
M yo (0) - H y3 (0) 
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" ^xy3(®) 


(11*7) 


M xyl (O) - M X y 2 (O) 

01x <0) “ ^2x(°> 

M y0 (O) - My 3 (0) (II*®) 

/?0y(O) ~ ^3y(0> 

M xy o(0) - M xy3 (0) 

^Ox(°) " ^3x(0> 

Enforcing equations (II-6) and (II-7) In the governing equations 
yields the following: 

f A 12 a 22 a 2s 1 f e 1 f"yil f"yu 1 


A 1 A 1 A 1 

A 12 A 22 A 26 


A 1 A 1 A 1 
16 26 66 


A° A ° A° 
A 12 A 22 26 


A° A ° A° 
16 A 26 66 


* N yU 

+ * N xylj s j G J “ 0 

N yU 

N xylj 


(II-9) 


fil 12 B 22 B 26 L 


A 1 B 26 B 66 


% n 


M yJj s j 

M xyij s j + B l s c a j 


(11-10) 


( j - 1- 4 ) 
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The expressions in (II-9) and (11-10) are defined below 


r ^ 

c y 


<■ 1 

Cv 

£ + 

c v 

Yx r 

k 4 


Cu 

Ik 4 


C u 

k 4 


* * * 2 
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Cv 

_ 1 

* 

a 26 

* 

A16 

* 

* a 66 

* 1 
a 12 

C u 

J 

A 

* 

a 26 

* 

a 12 
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*< 

• 

•k 

a 16 
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* * 

Cy- 

KM 

_ 1 

a 66 (Nyl + - A2$ (N X yl + N X yo) 

Cu 

k 4 

A 

-A 26 (Nyl + Nyo^ + a 22 (N X yl + N X yo) 

k * 


(II-H) 

(H-12) 


(H-13) 


( 11 - 14 ) 


Aj^ and A]_ are functions of Ay , By and Dy . The superscripts * 
implies a summation of the upper and lower sublaminates. 


Continuing with the derivation one can substitute the expressions set 
forth into equation (1-7) as well as 10 and 11 in the report. This gives the 
following expression for the total energy release rate. 


G 


1 

2 


b 

r 

f h/2 

■ 

I N x £« 

•4 

<h/2 


DY 


( 11 - 15 ) 


The concept of free -expans ion in the x-direction is implemented to 
find the strain induced by the non-mechanical effects on the structure. 
Setting N x — O for each ply in Equation (1-5 ) and using the boundary 
conditions of (II-2), (II-4) and (II-6) allows the following. 



- -1* / 

c c 


k 

c 

c 


( 11 - 16 ) 
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where 

T k - h k cf q £ 2 ♦ h k cT - (C) 

E* - h k (Qn + Cv Q12 + Cu Qi 6 ) 

T k " Ql2 fl” h k + Q k 6 F^ h k -t- B k 2 F3** - (n2 M ) k 

E k - Q k x h k + Q k 2 h k Cd n + Q k 6 h k Cd k 2 + B 12 Cd 31 (II -17) 

NM 

Superscript k represents the ply. Expressions Cdy and are found 

by substituting the conditions of (II-2) into (II-3). 

r 1 A i „i V 1 f .1 -i 

a 22 a 26 b 22 a 12 b 26 

[cdl- - A26 a 66 b 26 a 16 b 66 

1 J 3x2 

b 22 b 26 d 22 b 12 d 26 

1 1 1 T* f 

a 22 a 26 b 22 N y 

f”” ’ “ a 26 a 66 d 22 ‘ N xyl 

j x 1 111 

b 22 b 26 d 22 l M yl 

Sublaminate 3 has By - 0 due to symmetry of the structure. When 

J NM NM 

considering these plies, the term F 3 and F2 are found by substituting the 

boundary conditions (II-4) into (II-5). 
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This gives a second expression of for this sublaminate 




NM 
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•m 

t * 

2x1 - 


A26 


A26 


A66 


yl 

r i 


NyO 
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NM 


( 11 - 20 ) 


To find the total strain associated with the non-mechanical effects, 
it is necessary to sum the force over the entire structure and set it to zero. 
These are used in order to obtain Equation (16), (17) and (18) in the report 
on page 17. Substituting this in Equation (11-15) gives the total strain 
energy release rate expression per unit length 


G 


^.(^ T +6 (* T - 

2 K C 


*:> 


(B* 

u 


c T 4- I*) 
u 



) 


The expression c* - is 

c.u 

that ply. 


in essence the total mechanical strain of 


INTERLAMINAR STRESSES 

The interlaminar stresses of the structure are defined in Equations 
(8) and (9) 

fc x " ^xy,y ” ^xyij sj Gj e 

t y “ N y*y " N yU S J 2 G J e " sjy (j “ 3 ~ 4 ) (n-2i) 

While sj , the positive roots reulsting from the ploynominal 
Eg s 8 + Eg s$ + E 4 s 4 + E 2 s 2 + E 0 0, (11-22) 

are independent of the hygrothermal effects, the rest of the terms are not. 

Solving Equations (II-9) and (11-10) gives the term Gj while N X yjj an d 
Nylj are found from 
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'W 


A 22 

a 26 
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B 2 2 

b 26 



1 

1 

1 

1 

N xyJ 

k ■* 


A 26 

a 66 

b 26 

B 66 


(H-23) 


(j - 1,2, 3, 4) 


where vj , Uj and aj are found by imposing the boundary conditions on the mode 
shapes. They are dependent on the four values of sj as well as the 
sublaminate stiffness matrices. 
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PROGRAM START (INPUT, OUTPUT, TAPE5=«NPUT, TAPE 6 -OUTPUT) 

THIS PROGRAM IS FOR THE FINAL PAPER 8-16-87 
DIMENSION STATEMENTS 

REAL BG(C),E(9),GG(1|,C). MATR32(3,2), MATR33(3.3). MATR 3 ( 3 ) , 

C STRAIN (25), SAVE3 (3) » SAVE 33 (3, 3), WKAREA (99) , ZR 
COMPLEX SJ ( 8 ) 

DOUBLE PRECISION BNEG, A, C, DIFF1, DIFF2, UNSY (2) , UNSX (2) , 

C SSSS, SSSC, SSCC, SCCC, ZZZ0(0:50) ,J22,J26,FNM(3) .S1NM.S2NM, 

C MEMSY, MEMSX, F11M, F22M, SSI, SS 2 , SSY, SSX, ZZZ1(0:50), 

C TH I CK (1+0) .THETA (50) , El (50) . E2 (50) , CCCC, HSS (5) ,HSN1 ,HSN 2 , 
c Q ( 6 , 6 , 50 ) , ZO (0:40) , A 0 ( 6 , 6 ), A 1 ( 6 , 6 ) , U 13 , UU, BG 1 , BG 2 , 

C NXO(C), E15.E16.E17, E 18, E19, ZTT(0:C0), 

C NY 1 (C) , X (2) , Y (3) . CV, CU, W(2), ZZ (3) , VV 11 , VV12 
OOUBLE PRECISION ALPHA (C) , PHI (4), GAMMA (C) , NX1 (A) , 

C B 1 ( 6 , 6 ) , BO ( 6 , 6 ) , DO ( 6 , 6 ) , D1 ( 6 , 6 ) , F(C,C), VV13, WlC, J 66 , 

C RDLT , RTA1, RTA2 , RSB1, RSB2 , U12, Ull, A1NM, 

C NXY1(C), MY1(C), MXY1 (U) , WD (2,3) , CD (3, 2), WIDTH, 

C V12(50). V2 1 (50) , SS, CC, K 66 , K26, K22, 

C Z1 (0: AO) , FX, FY, G I (2,35) . K16.K12, H 66 

DOUBLE PRECISION SV(5), SU (5) . AL, SC, S ( 8 ) , DY, 

C G 1 2 (50) , G31 (50) , C2 , Cl. THETV, THETU, Gill (2,35). CS, 

C DEL, HO, HI, H22, HE, HG, HNY (50) , HNXY (50) , HM3.CVNM, 

C C11, C12, C22, C26, CCC, D, C55. C 66 , H 26 , SQ, DUM.CUNM, 

C CONY, CONXY, SMNY, SMNXY, 

C SB1 ,SB2,TA1 ,TA2,ATHMl (CO) ,ATHM2(C0) ,ATHM6(40) ,BSW2(C0) , 

C NMNYO , NMNXYO , G II (2 , 35) . 

C DVV11, DVV12, 0U13. OUlC, DF (C.C) , DX, ATH , CCONY , CCONXY , 

C BSW 6 (1*0) , DELTEMP, BSW1 (1*0) , CMOIST(35). 

C SIGX (0:A0, 79:120) , 

C NMNY1 , NMNXY 1 , NMMX 1 , NMMY 1 , S I G Y (0 : CO , 79 : 1 20) 

DOUBLE PRECISION NMSTO ( 50 ) , NMST2 (50) , NMST 3 ( 50 ) , TNC.UNCL, 

C T1 (50) , T12(50), T1 3 (50) , EX (50) , EXX (50) , EX3 (50) , JY, 

C EXNC, ESTAR, TSTAR, TNMST, GLC(0:50), NXNM(50) ,B12 (50) 

DATA Q/l800*0.0/, ZR/O.OO/ 

************************************************£************************ 

A************************************************************************ 

DATE OF PROGRAM : SEPT. 1, 1 987 

THIS IS THE FINAL PROGRAM FOR PREDICTING THE ENERGY RELEASE RATE OF 

COMPOSITE LAMINATES INCLUDING HYGROTHERMAL EFFECTS 

HOWEVER, IT ONLY CONSIDERS EXTENSION EFFECTS OF STRAIN WHEN DEALING 

WITH HYGRALTHERMAL EFFECTS LIKE WHITTNEY'S PAPER 

EXCEPT ON A PLY BY PLY ANALYSIS BASIS OF THE HYGRALTHERMAL EFF-ECTS 


THE INPUT ALLOWS FOR: THE LAMINATE LAY-UP TO CHANGE AND POSITION 
OF THE CRACK, DIFFERENT STRAIN VALUES TO BE EVALUATED, 
(UP TO AO LAMINATES AND 25 DIFFERENT STRAIN VALUES) 

AND FOR THE EVALUATION OF ONE MOISTURE CONSTANT OR A 
RANGE OF THE MOISTURE CONSTANT FROM 0 TO 1.2. 

THIS PROGRAM IS FOR THE GIVEN DATA TO BE IN ENGLISH UNITS. 

ALL LAMINATES ARE EVALUATED WITHOUT THRMAL EFFECTS AUTOMATICALLY 


************ ***************** ******************************************* 
*************************************** 3 ********************************* 


LZQ IS THE NUMBER OF DIFFERENT LAMINATE (OR CRACK POSITIONS) 
TO BE EVALUATED. 


RE AO (5,*) LZQ 

DO COO L ZZ = l.LZQ 


ORIGINAL PAGE S3 

OF POOR QUALITY 



RE AO (5 . *) WIDTH, NPLYO, NPLY1 , AL 
NEXTPL = NPLYO + 1. 

TPLY = NPLYO + NPLY1 

FOR EACH PLY IN THE SUBLAMINATE, THE MATERIAL CHARACTERISTICS 
MUST BE READ IN. 

PI = 1*. * ATAN(1.) 

HO - 0.0 
ATH *=0.0 
H 1=0.0 

DO 3 LK = 1, TPLY 
20 (LK) *= 0.0 
3 Z1 (LK) = 0.0 


DO 5 I - 1, NPLYO 

READ (5.*) THICK(I), THETA (I ) , El(l), E2(l) 

READ (5,*) VI 2(1) 

READ (5.*) G12(l), G3U0 
THETA ( I ) = THETA (I) * PI / 180. 

5 HO = HO + THICK (I ) 

DO 10 I « NEXTPL, TPLY 

READ (5.*) THICK(I), THETA (I) , E1(l), E2 (I) 

READ (5,*) VI 2 (I ) 

READ (5»*) G12 (I) , G3 1 (I) 

THETA (I) = THETA (I) * PI / 180. 

10 HI - HI + THICK(I) 

********************************************************************* 
********************************************************************* 
THESE ARE WRITE STATEMENTS TO CHECK THE INITIAL CONDITIONS OF THE 
SUBLAMINATE AND VALUES READ IN 

********************************************************************* 

********************************************************************* 

EACH PLY MAY HAVE DIFFERNT PROPERTIES SO THE PROPERTY OF EACH 

WRITE ( 6 , 289 ) 

CC 

WRITE (6,201) WIDTH 
WRITE (6,202) NPLY1, NPLYO 

WRITE (6,20*4) 

DO 15 J - 1, TPLY 
JJ = TPLY + 1 - J 
WRITE (6,206) J 

WRITE (6,207) THICK(JJ), THETA (J J) *1 80/P I 
WRITE (6,208) El (JJ) /1E+06 , E2(JJ)/lE+06 
WRITE (6,209) VI 2 (J J) 

15 WRITE (6,2114) G12 (JJ) /I E+06 , G31 (JJ) /1E+06 

******************************************************************** 
******************************************************************** 
DETERMINE THE 2 COMPONENT OF ALL LAMINATES 

CHECK = 0.00000001 
ZTT(O) =0.0 
Z0(0) = -HO / 2.0 
DO 20 I = 1, NPLYO 

ZTT(I) = THICK (I) + ZTT(l-l) 

20 20(1) = THtCK(l) + ZO(l-l) 


Z1 (NPLYO) = -HI / 2.0 
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DO 25 I = NEXTPL ,TPLY 

ZTT(I) - THICK(I) + ZTT(l-l) 

Z1 (I) - THICK(I) + Z1 (1-1) 

********************************************************************* 
FIRST READ IN THE NUMBER OF STRAINS TO BE EVALUATED AND THEIR VALUE 
THEN READ IN IF THE MOISTURE CONTENT SHOULD VARY OVER 0 TO 1.2 OR 
BE A CONSTANT. 

NSTRA = NUMBER OF VARIOUS STRAIN VALUES 

IF MOISTV - 1 ... CMOIST VARIES OVER 0 TO 1.2 
IF MOISTV = 0 ... CMO I ST IS A SPECIFIC VALUE 
********************************************************************* 
READ (5,*) NSTRA 
DO 27 J=l, NSTRA 
27 READ (5,*) STRAIN(J) 

DO 400 1ST - 1, NSTRA 
READ (5,*) RDLT.RSBl .RSB2.RTA1 ,RTA2 

WRITE (6,231) STRAIN (LST) ,RDLT, RSB1, RSB2, RTA1, RTA2 

READ (5,*) MOISTV 
IF (MOISTV. EQ.O) READ (5,*) CM 
IF (MOISTV. EQ.O) MMC - 1 
IF (MOISTV. EQ.O) WRITE (6,232) CM 
IF (MOISTV. EQ.1) MMC - 25 
IF (MOISTV. EQ.1) WRITE (6,233) 

********************************************************************* 
DO 300 JM - 1 ,MMC +1 

FIND Q'S AS WELL AS Q-BAR , SAVING Q-BAR 
AND READ AND CALCULATE THE HYGRO THERMO EFFECTS 

********************************************************************* 
DO 200 IZZ - 1,2 
III - 0 

IF (IZZ. EQ.1) JMM - JM 
IF (I ZZ.EQ.2) JMM - 0 
I F (JM.EQ. 1 . AND .IZZ. EQ.1) LIL - 1 
IF (JM.GT.l .AND. IZZ.EQ.2) GO TO 200 

NMNY1 - ZR 
NMNXY1 = ZR 
NMMX1 - ZR 
NMMY1 - ZR 

NMNYO - ZR 
NMNXYO « ZR 
HM3 - ZR 
SMNY - ZR 
SMNXY = ZR 

DO 24 1=1,5 
E(l) = ZR 
E (1+4) = ZR 
24 HSS(I) = ZR 

DO 26 1=1,6 
DO 26 J= 1,6 
NXNM ( I ) = ZR 
DF (I , J) = ZR 
AO(l,J) = ZR 
BO (I , J) = ZR 
DO ( I , J) = ZR 
A1 (I ,J) = ZR 
B1 (I , J) = ZR 
D1 (I ,J) = ZR 
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DO 28 MM * 1 ,TPLY 

28 I F ( THICK (MM) .GT. ATH ) ATH = THICK (MM) 

IZ2 = 2 IS FOR LAMINATE WITHOUT ANY HYGROTHERMAL EFFECTS 
IZZ * 1 IS FOR HYGROTHERMAL EFFECTS CONSIDERED 

DO 30 1=1, TPLY 

READ THE HYGROTHERMO EFFECTS, BOTTOM PLY IS FIRST AND UPWARD 

IF (I ZZ.EQ.2) GO TO 35 

IF (MOISTV.EQ.O) CMOIST(JM) = CM 
IF (M0ISTV.EQ.1) CMOIST(JM) = 0.05 * (JM-1) 
OELTEMP = RDLT 
SB1 = RSB1 
SB2 = RSB2 
TA1 = RTA 1 
TA2 ■= RTA2 
GO TO AO 

35 DELTEMP ■= ZR 

CMOIST(JM) = ZR 
SB1 ■= ZR 
SB2 = ZR 
TA1 = ZR 
TA2 = ZR 

AO V21 (I) - V12 (I) * E2 ( I ) / El (I) 

Cl 1 = El (I) / (1 - V12 (I) * V21 (I) ) 

C 1 2 = E2 ( I ) * VI 2 ( I ) / ( 1 - VI 2 (I ) * V2 1 ( I ) ) 

C22 = E2(l) / ( 1 - VI 2 ( I ) * V21 (I) ) 

CAA = G31 (I) 

C55 - G31 (I) 

C66 = G 1 2 ( I ) 

SS = DSIN (THETA (I )) * DS I N (THETA ( I ) ) 

CC = 1 - SS 

CS - 0.5 * DS I N (2*THETA (I ) ) 


SSSS = SS * SS 

SSSC = SS * CS 

SSCC = SS * CC 

SCCC = CC * CS 

CCCC = CC * CC 

Qd.l.l) - Cll * CCCC + 2 * (Cl 2 + 2 * C66 ) * SSCC 
c + C22 * SSSS 

QO.2,1) « (Cll + C22 - A * C66) * SSCC + C12 * ( SSSS 
C + CCCC ) 

Q (2, 2 , I ) = Cll * SSSS + 2 * ( C12 + 2 * C66) * SSCC 
C + C22 * CCCC 

Q(l,6,l) - (Cll - C12 - 2 * C66) * SCCC 
C + (C 1 2 - C22 + 2 * C66 ) * SSSC 

Q(2,6, I) = (Cll - C12 - 2 * C66 ) * SSSC 
C + (C12 - C22 + 2 * C66 ) * SCCC 

Q (6,6, I) « (Cll + C22 - 2 * C 1 2 - 2 * C66) * SSCC 
C + C66 * ( SSSS + CCCC ) 

Q(A,A, l)= CAA * CC+ C55 * SS 
Q (5.5. I) = CAA * SS + C55 * CC 
Q(A,5,I) = CS * (CAA - C55) 

Q (6.2. 1) = Q (2.6. I) 

Q(6,l,l) = Q (l ,6, l) 

Q (2, 1 , 1) = Q (1 ,2, I) 


HSS(l) = HSS(l) + Q (1 , 2, I ) 
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HSS (2) - HSS (2) + Q (2, 2 , I ) 
HSS (3) - HSS (3) + Q (2,6, 1) 
HSS (i*) « HSS(4) + Q (1 ,6, 1) 
HSS (5) « HSS (5) + Q (6.6, I ) 


ATHM1 (I) - TA1 * CC + TA2 * SS 
ATHM2 ( I ) = TA1 * SS + TA2 * CC 
ATHM6(I) * CS * ( TA2 - TA1) 

BSW1 (I) *= SB 1 * CC + SB2 * SS 

BSW2 ( I ) * SB1 * SS + SB2 * CC 

BSW6(I) « CS * ( SB2 - SB1 ) 

30 CONTINUE 

FIND THE A, B, AND D MATRICES FOR THE LOWER AND UPPER SUBLAMINATE. 
ALSO FINDS HYGRALTHERMAL EXPRESSIONS ON A PER LAMINA ANO PER 
SUBLAMINATE BASIS. 


ZZZO (0) = Z0(0) * Z0(0) * ZO (0) 
DO 45 I - 1, NPLYO 


NXNM(I) - (Q(l,1,l)*( ATHM1 (I) *DELTEMP + BSW1 (I) *CMO JST(JM) ) 

C + Q (1 ,2, I) * ( ATHM2 (I) *DELTEMP + BSW2 ( I ) *CMO I ST (JM) ) 

C + Q(l,6.l)*( ATHM6(l)*DELTEMP + BSW6 ( I ) *CMO I ST (JM) ) ) * THICK(I) 


NMNYO= NMNY0+(Q(1,2, l)*( ATHM1 ( I ) *DELTEMP + BSW1 (l ) *CMO I ST (JM) ) 

C + Q (2 , 2 , I) * ( ATHK2 ( I ) *DELTEMP + BSW2 (I) *CMOI ST (JM) ) 

C + Q(2.6, l)*( ATHM6 ( I ) *DELTEMP + BSW6 ( I ) *CMO 1 ST (JM) ) ) * THICK (I) 

NMNXYO* NMNXY0+(Q(1,6, l)*( ATHM1 (I ) *DELTEMP + BSW1 (I) *CMOIST (JM) ) 
C + Q (2,6, I) * ( ATHM2 ( I ) *DELTEMP + BSW2 ( I ) *CM0 1 ST (JM) ) 

C + Q (6,6, I) * ( ATHM6(I)*DELTEMP + BSW6 ( I ) *CM0 1 ST (JM) ) ) * THICK(I) 


ZZZO(I) = Z0(l) * ZO (I ) * Z0(l) 

B 1 2 (I ) = Q ( 1 , 2 , I ) *0 . 5* ( (ZO ( I ) *Z0 ( I ) ) - (ZO ( I - 1 ) *Z0 ( I - 1 ) ) ) 

DO 45 L -1,6 
DO 45 J *= 1,6 

I F ( REAL ( Q (J, L , I ) ).EQ.ZR) GO TO 45 
AO (J , L) - AO ( J , L) + Q (J , L, I ) * THICK(I) 

BO(J.L) - BO (J , L)+Q (J , L, I) *0.5* ( (ZO (I) *Z0 (I) ) - (ZO (I -1) *Z0 (I -1) ) ) 
DO (J , L) - D0(J,L)+Q(J,L. l)/3-0*( ZZZO(I) - ZZZO(i-l) ) 

45 CONTINUE 

ZZZ1 (NPLYO) « ZI (NPLYO) * Z1 (NPLYO) * Z1 (NPLYO) 


DO 50 I * NEXTPL.TPLY 

ZZZ1 (I) = Z1 (I) * Z1 (I) * Z1 (I) 


NXNM(I) = ( Q(1, !,!)*( ATHM1 ( I ) *DELTEMP + BSW1 ( I ) *CMO I $T (JM) ) 

C + Q(l,2,l)*( ATHM2 ( I ) *DELTEMP + BSW2 ( I ) *CMO I ST (JM) ) 

C + Q (1 ,6, 1 ) * ( ATHM6 (I ) *OELTEMP + BSW6 ( I ) *CMO I ST (JM) ))* THICK(l) 


NMNY1= NMNY1+ ( Q(l,2,l)*{ ATHM1 ( I ) *DELTEMP + BSW1 ( I ) *CMO I ST (JM) ) 
C + Q (2,2, I) * ( ATHM2 (I) *DELTEMP + BSW2 (I) *CMO I ST (JM) ) 

C + Q (2,6, I) * ( ATHM6(I)*0ELTEMP + BSW6 ( I ) *CMO I ST (JM) ) )* THICK(I) 
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NMNXY1- NMNXY1+(Q(1,6,I)*( ATHM1 (I ) *0ELTEMP + BSW1 (I ) *CMO I ST (JM) ) 
.C + Q(2,6,l)*( ATHM2 (I) *DELTEMP + BSW2 (I ) *CM0 1 ST (JM) ) 

C + Q(6,6,I)M ATHM6(I)*DELTEMP + BSW6 (I ) *CMO I ST (JM) ) )* THICK(I) 


NMMX1- NMMX1 + 0.5 * ( Z 1 ( I ) *Z 1 ( I ) - Z1 (I -1) *Z1 (I -1) ) * 

C ( Q (1.1, DM ATHM1 (I) *DELTEMP + BSW1 (I) *CM0I ST (JM) ) 

C + Q(1.2,l)*( ATHM2 (l)*DELTEMP + BSW2 (I ) *CM0 1 ST (JM) ) 

C + Q (1 ,6, I) * ( ATHM6(I)*DELTEMP + BSW6(I)*CM0IST(JM) ) ) 


NMMY 1 = NMMY1 + 0.5 * ( Z1 (l)*Z1 (I) - Z1 (I-1)*Z1 (1-1) ) * 

C ( Q(l,2,l)*( ATHM1 (I) *DELTEMP + BSW 1 ( I ) *CM0 1 ST ( JM) ) 

C + Q(2, 2, l)*( ATHM2 (I ) *DELTEMP + BSW2 (I ) *CM0l ST (JM) ) 

C + Q(2,6, DM ATHM6(I)*DELTEMP + BSW6 ( I ) *CM0 1 ST (JM) ) ) 


B 1 2 ( I ) - Q ( 1 , 2 , I ) *0 . 5* ( (Z 1 ( I ) *Z 1 ( I ) ) - (Z 1 ( I - 1 ) *Z 1 ( 1 - 1 ) ) ) 

DO 50 L«1,6 

DO 50 J«1,6 

IF ( REAl( Q(J,L, I) ).EQ.O ) GO TO 50 
A 1 (J, L) -A1(J,L) + Q(J,L,I) * THICK(I) 

B 1 (J , L) - B1 (J,L)+Q(J,L, 0*0.5* ((Zl (I)*Z1 (I))-(Z1 (1-1) *Z1 (1-1))) 
D 1 (J , L) - 01 (J,L)+Q(J,L,l)/3.0*( ZZZ1 (I) - ZZZl(l-l) ) 

50 CONTINUE 


****************** ************************************************** 


SEE IF COUPLING IS TAKING PLACE 

******************************************************** *********** 
COUPL = 2 
00 60 1-1.6 
DO 60 J*l,6 

IF ( REAL (BO (I ,J)) .GT. CHECK ) COUPL-1 
60 IF ( REAL (B 1 (I ,J)) .GT. CHECK ) COUPL-l 

IF (REAL (D 1 (2,6) ) .GT. CHECK) COUPL-1 
IF ( REAL (00(2,6)) .GT. CHECK ) COUPL-1 

IF ( COUPL. EQ.l .AND. LIL.EQ.1 ) WRITE (6,205) 

IF ( COUPL. EQ. 2 .AND. LIL.EQ.1 ) WR!TE(6,210) 
****************************************************************** 


CHECK THE SIGN OF THE PEEL STRESS *************************** 

HSN 1 - NMNY1 + NMNYO 
HSN2 = NMNXY1 + NMNXYO 

****************************************************************** 
HDD- HSS (2) * HSS(5) - HSS (3) * HSS (3) 

HE - HSS (3) * HSS (A) - HSS (1) * HSS (5) 

HE - HE /HDD 

HG - HSS (1) * HSS (3) - HSS (2) * HSS (A) 

HG - HG / HDD 
DO 65 l-l.TPLY 

HNY(I) = ATH * STRAIN(LST) * ( Q(l,2,l) + Q(2,2,l) * HE + 

C Q (2 ,6, I) * HG ) 

HNXY(I) = ATH * STRAIN(LST) * ( Q(l,6,l) + Q(2,6,l) * HE + 
C Q (6 , 6 , I ) * HG) 

65 CONTINUE 


DO 70 I = l.NPLYl 


UUUUNUU 


H«3 * HM3 + ATH * HNY(l) * ( NPLY1 - I + .5 ) 
SMNY « SMNY + HNY(I) 

70 SMNXY «= SMNXY + HNXY(I) 


IF ( HM3.GT.ZR) GO TO 85 

IF (LIL.EQ.1) WRITE(6,*)' CASE OF COMPRESSIVE PEEL STRESS 
WRITE (6,218) 

DO 75 l-l.TPLY 

WRITE (6,220)THETA(I) ,HNY(I) ,HNXY(I) 

WRITE (6,*) ' THE MOMENT CALCULATED WAS = ' ,HM3 

85 DO 80 1*1,6 
DO 80 J*1 ,6 

IF( ABS ( REAL (BO (I , J) ) ) .LT. CHECK) BO(l,J)«=ZR 
80 IF ( ABS ( REAL (B1 (I ,J)) ) .LT. CHECK) B1(I,J) -ZR 

******************************************************************** 
********************************************************* 

DEFINE SOME PARAMETERS NEEDED IN THE PROGRAM 

********************************************************* 

H22 - B1 (2.2) + HI / 2.00 * A1 (2,2) 

H26 - B1 (2,6) + HI / 2.00 * A1 (2,6) 

H66 - B1 (6,6) + HI / 2.00 * A1 (6,6) 

C22 - BO (2 , 2) + HO / 2.00 * A1(2,2) 

C26 - BO (2,6) + HO / 2.00 * A1 (2,6) 

C66 - 80(6,6) + HO / 2.00 * A1 (6,6) 

K22 * A 1 (2,2) + AO (2,2) 

K26 - A1 (2,6) + AO (2,6) 

K66 ■ A1 (6,6) + AO (6,6) 

K12 - A1 (1,2) + AO (1 , 2) 

K16 - A1 (1,6) + AO (1 , 6) 

D - K22 * K66 -( K26 * K26) 

E 15 - DO (2,2) - H0/2*B0 (2 , 2) 

E 1 6 - DO (2,6) - H0/2*B0 (2 , 6) 

£17 - DO (6,6) - H0/2*B0 (6,6) 

E 1 8 - BO (1,2) - H0/2*A0 (1,2) 

E19 - BO (1,6) - H0/2*A0 (1,6) 

VV11 - ( K26 * H26 - K66 * H22 ) / D + ( HI / 2.00 ) 

VV12 * ( K26 * C26 - K66 * C22 ) / D + ( HO / 2.00 ) 

VV13 - ( K26 * H66 - K66 * H 26 ) / 0 

VVU 4 « ( K26 * C66 - K66 * C26 ) / D 

U11 ■ ( K26 * H22 - K22 * H26 ) / D 

U1 2 *= ( K26 * C22 - K22 * C26 ) / D 

U13 *= ( K26 * H26 - K22 * H66 ) / D + ( HI / 2.00 ) 

U14 *= ( K26 * C26 - K22 * C66 ) / D + ( HO / 2.00 ) 

F (1,1) - 01(2,2) + B 1 (2,2) * HI / 2.0 + H22*VV11 + H26 * Ull 
F (2,1) * H22 * VV1 2 + H26 * U12 

F (3, 1) - D 1 (2,6) + B 1 (2,6) * HI / 2.0 + H22*VV13 + H26 * U1 3 
F (A, 1) - H22 * VV14 + H26 * U14 

F (2,2) = DO (2 , 2) - BO (2,2) * HO / 2.0 + C22 * VV12 + C26*U12 
F (3,2) = H26 * VV1 2 + H 66 * U12 

F (1+ , 3 ) = H26 * VV1L + H 66 * UU 

F (3,3) - D1 (6,6) + B 1 (6,6) * HI / 2.0 + H26*VV13 + H66 * U 1 3 

• F (4 , 2) = 00(2,6) - BO (2,6) * HO / 2.0 + C22*VV14 + C 26 * U14 
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c - 3 - 



ooo o ooo o ooo 


F(A,A) - DO ( 6 , 6 ) - BO( 6 , 6 ) * HO / 2.0 + C 26 ftVVlA + C 66 ft uiA 
DX - K22 * K 66 

DVV 11 - - K 66 * H22 / DX + ( HI / 2.00 ) 

DVV12 *= - K 66 * C22 / DX + ( HO / 2.00 ) 

DU13 « - K22 * H 66 /OX + ( HI / 2.00 ) 

DU 1 i* « - K22 * C 66 / DX + ( HO / 2.00 ) 

DF (1,1) *= D 1 (2,2) + B 1 (2,2) * HI / 2.0 + H22*DVV1 1 
DF (2,1) « H22 * DVV12 

DF (2,2) - DO (2 , 2) - 80(2,2) * HO / 2.0 + C22 * DVV12 
DF (A, 3) - H 66 * DU1I* 

OF(3,3) « Dl( 6 , 6 ) + B1 ( 6 , 6 ) A HI / 2.0 + H 66 * DU13 
DF (I* , i*) «= DO ( 6 , 6 ) - 80(6,6) * HO / 2.0 + C 66 * Dim 


W(1) « F(3,3)*( F (2,2) *F (A , I*) — F (A, 2 ) *F (4,2) )- F ( 3 , 2 ) *f ( 3 . 2 ) * 
C F(A,A) + 2*F(i»,3)*F(ii.2)*F(3.2) - F (2,2) *F (1* , 3) *F (fc, 3) 

W(2) « - F (3,3)*(F (2,2)*A0(5.5) + F (A,A) *A0 (I* , 1») - 2*F (I*. 2) 

* AO(l»,5) ) ~ A1 (5,5) * ( F (2,2) * F(A,A) - 
F (A, 2 ) * F (A, 2) ) + F(3,2) * F( 3 , 2 ) * AO (5.5) * 

2.0 *F(A,3)*F(3,2)*AO(A,5) + F(A,3) * F(A,3) * AO (A, A) 


X (1) -F (3,1) * (F (2,2) *F (A, A) -F (A, 2 ) *F (A, 2 ) ) -F (3,2) * (F ( 2 , 1 ) *F (A, A) 
~ F (A , 1 ) *F (A , 2) ) + F (A,3) * (F (2,1) *F (A, 2 ) - F (A, 1 ) *F (2, 2) ) 

X (2) « - F (3, 1) * ( F (2,2) * AO (5,5) + F (A, A) * AO (A, A) 

- 2 * F (A, 2 ) * AO (A, 5) ) - AT (A, 5) * ( F ( 2 , 2 ) * F(A.A) 

- F (A, 2) * F (A, 2) ) + F (3,2) * ( AO (5,5) * F (2,1) - 
F (A, 1) *A0 (A,5) ) - F (A , 3) * ( F (2,1)*AO(A,5) - F (A, 1) *A0 (A, A) ) 


Y(l) * F (3, 1 ) * (F ( 3 , 2 ) *F (A , A) - F(A, 3 )*F(A, 2 ))- F (3,3) * (F ( 2 , 1 ) * 
F (A, A) - F(A,1)*F(A,2))+F(A,3)*(F(2,1)*F(A,3)-F (A,1)*F (3,2)) 

Y (2) - 0 - A 1 (A,5)*(F ( 3 . 2 ) *F (A, A) - F (A, 2 ) *F (A, 3)) 

“ F(3,l) * ( F (3,2) *A 0 (5.5) " F(A,3) * AO (A, 5) ) + 

A1 (5.5) * ( F ( 2 , 1 ) * F (A, A) - F (A , 1) * F (A, 2) ) + 

F (3,3) * ( F (2 , 1 ) * AO (5,5) - F (A , 1) * AO (A, 5) ) 

Y(3) - A 1 (A,5) * (F (3,2) *A 0 (5,5) - F(A,3)*AO(A,5)) - A 1 (5,5) * 
c ( F ( 2 , 1 ) * AO (5,5) - F (A, 1 ) ft AO (A, 5 ) ) 


ZZ( 1 ) - F (3, 1) * (F ( 3 , 2 ) * F (A, 2 ) - F(A,3)*F(2,2)) - F (3.3) * (F (2, 1) * 

C F (A, 2) -F (A, 1 )*F (2,2)) +F ( 3 , 2 ) ft (F ( 2 , 1 ) *F (A,3)-F (A,1)*F (3,2)) 

ZZ ( 2 ) - F (3 , 1 ) * (F (A , 3) *A 0 (A , A) - F (3.2) *A0 (A,5) ) - A 1 (A , 5 ) * 

C (F (3.2) *F (A, 2 ) - F(A,3)*F(2.2))+ A1(5,5)*(F(2,D* F(A,2)- 

C F (A,1)*F ( 2 , 2 )) - F (3,3) * (F (A , 1 ) *A 0 (A , A) - F (2, 1) *A0 (A,5) ) 

ZZ (3) * 0 - A 1 (A , 5) * (F (A , 3 ) *A 0 (A , A) -F ( 3 , 2 ) *A 0 (A , 5 ) ) +A 1 (5,5) * 

C ( F (A, 1 ) ft AO (A, A) - F ( 2 , 1 ) * AO (A, 5) ) 

************************************************************************ 


NOW OBTAIN THE VALUES OF E SO THAT THE 8TH ORDER POLYNOMIAL MAY BE SOLVED 


C 

C 

C 


E (1) *= F (1,1) *W ( 1 ) - F (3.1)*X(1) + F (2, 1) *Y (1) - F(A,1)*ZZ(1) 

£ (3) = F ( 1 , 1) *W (2) - A 1 (A, A) ftW ( 1 ) - F (3, 1) *X (2) + A1(A,5)*X(1) 
+ F (2, 1) * Y (2) - F (A. 1) * ZZ (2) 

E(5) = (AO (A , A) *A0 (5,5) - AO(A,5)ftAO(A,5))*(F (1, 1)*F (3,3) 

- F ( 3 , 1 )*F ( 3 , 1 )) + (F (2 , 2 ) *A 0 (5,5) + F (A, A) *A0 (A, A) - 
2*F (A.2)*A0(A,5))ft(F (), 1)*A1 (5,5) - F (3. D *A 1 (A, 5) ) 



c - AT (It.lf) *W(2) + A1 (J.,5)*X(2)+ F (2,1)*Y(3) - F (4, 1) *ZZ (3) 

E(7) - - (AO (4,4)* AO (5.5) - AO (4.5) *AO (4.5) ) * (F (1 . 1) *A1 (5.5) 

C + F(3,3)*A1 (*.*)- 2*F (3. 1)*A1 (4,5)) " (A1 (4,4)*A1 (5.5)- 

C A1 (4,5) *A1 (4,5) ) * (F (2,2) * AO (5.5) + F (4, 4) * AO (4,4) 

C - 2 * AO (4,5) * F (4,2) ) 

E(9) = (AO (4,4) *AO (5,5) - AO(4.5)*AO(4,5) ) * 

C (A1 (4,4) * A1 (5,5) * A1 (4,5) * A1 (4,5) ) 

CALL UP SUBROUTINE TO SOLVE 8TH ORDER POLYNOMIAL 

NDEG - 8 
IER «= 0 

CALL ZPOLR (E.NOEG.SJ, IER) 

KK « 0 

******************************************************************** 

******************************************************************** 



IF 

(LIL.EQ. 

1) 

WR 1 TE (6 

,217) 



DO 

90 L * 

1. 

8 




S(L) - REAL 

(SJ (D) 




IF 

(REAL (SJ 

(L) 

|) .GT.O) 

KK - KK 

+ 1 

90 

IF 

(REAL (SJ 

(L) 

i) .GT.O) 

S (KK) - 

S(L) 


DO 95 kk 

SE 

1,4 



95 


IF (LIL. 

EQ 

.1) WRITE (6,221) 

KK, S (KK) 


******************************************************************** 


******************************************************************** 

NOW FIND THE UNCOUPLED S VALUES AND THOSE OF THE MEMBRANE 

BNEG- OF (1,1) * AO (4,4) + OF (2,2) * A1 (4,4) 

A « OF (1,1) * OF (2,2) - DF (2,1) * DF (2,1) 

C - AO (4,4) * A1 (4,4) 

SQ * DSQRT (BNEG * BNEG - ( 4.0 * A * C) ) 

DIFF1 - DABS ( BNEG - SQ ) 

DIFF2 = DABS ( BNEG + SQ ) 

IF (DIFF1 .GT.DIFF2) GO TO 100 

UNSY(l) = DSQRT ( (BNEG+SQ) / 2.0 / A ) 

GO TO 105 

100 UNSY(l) - DSQRT ( (BNEG-SQ) / 2.0 / A ) 

105 UNSY (2) = DSQRT ( (BNEG/A) - UNSY(l)* UNSY(1) ) 

BNEG* DF (3, 3) * AO (5.5) + DF(4,4) * A1 (5,5) 

A - DF (3.3) * DF (4,4) - DF(4,3) * DF(4,3) 

C - AO (5,5) * A 1 (5,5) 

SQ - DSQRT (BNEG * BNEG - ( 4.0 * A * C) ) 

0IFF1 « DABS ( BNEG - SQ ) 

DIFF2 « DABS ( BNEG + SQ ) 

IF (DIFF1.GT.DIFF2) GO TO 110 
UNSX(I) * DSQRT ( (BNEG+SQ) / 2.0 / A ) 

GO TO 115 

110 UNSX(l) = DSQRT ( (BNEG-SQ) / 2.0 / A ) 

115 UNSX (2) - DSQRT ( (BNEG/A) - UNSX(l) * UNSX (1) ) 

******************************************************************** 
IF (LIL.EQ.1) WRITE (6,224) UNSX(l), UNSX (2) 

IF (LIL.EQ.1) WRITE (6,223) UNSY(l), UNSY (2) 

******************************************************************** 


NOW THE S FOR THE MEMBRANE ONLY 



c 


71 


F11« • HI/2. * ( AT (2,2) * (0VV1 1 + DVV1 2*HO/H 1*A 1 (U , U) /AO (A ,U) ) 

+ B 1 (2,2) ) 

F22M - Hl/2. * ( A 1 (6,6) * (U 1 3 + UU»*H0/H1*A1 (5.5) /AO (5,5) ) 

C + B1 (6,6) ) 

MEMSY = DSQRT ( Al(l4,l*) / F 1 1M ) 

MEMSX «= DSQRT ( A1 (5,5) / F22M ) 

IF (LIL.EQ.1) WRITE (6,219) MEMSX, MEMSY 

******ft*******^******rt*****rtrt***>V*rt**ftftrt*rtrt*}V:fr*rtrt*ft*********rtrt*j'c*rt}'; 

OUM - AO (2,2) * AO (6,6) - AO (2,6) * AO (2,6) 

S1NM - ( AO(6,6)*NMNYO - AO (2,6) *NMNXYO ) / DUM 
S2NM - ( AO (2,2) *NMNXYO - AO (2,6) *NMNYO ) / DUM 
IF (LIL.EQ.1) WRITE(6,288) S1NM.S2NM 

IF ( COUPL.EQ.2) GO TO 130 

SOLVE FOR WD, CD, CU AND CV 

WD(l.l) - (AO (2 , 6) * AO (1 ,6) - AO (6,6) * AO 0.2) ) / DUM 

WD (1 , 2) - (AO (2,6) * BO (2,6) - AO (6,6) *B0 (2,2) ) / DUM 

WD(1,3) - (AO (2,6) *B0 (6 , 6) - AO (6,6) *B0 (2,6) ) / DUM 
WD (2, 1) * (AO (2,6) * AO (1 , 2) - A0(2,2) * AO ( 1 , 6) ) /DUM 

WD (2,2) - (AO (2,6) *B0 (2,2) - AO (2,2) * BO (2,6) ) / DUM 

WD (2,3) - (AO (2,6) * BO (2,6) - AO (2 , 2) *B0 (6 , 6) ) / DUM 


MATR33 (1,1) - A1 (2,2) 
MATR33 (1,2) ■= A1 (2,6) 
MATR33 0.3) - B1 (2,2) 

MATR33 (2,1) «= MATR33 (1,2) 
HATR33(2,2) = A1 (6,6) 
MATR33(2,3) - B1 (2,6) 

MATR33 (3.1) - MATR33 (1,3) 
MATR33 (3,2) - MATR33 (2,3) 
MATR33(3.3) * D 1(2,2) 


MATR32 (1,1) - - A1(1, 2) 

MATR32 (1,2) = - B1 (2,6) 

MATR32 (2,1) = - A1 (1,6) 

MATR32 (2 , 2) •= - B1 (6,6) 

MATR32 (3,1) = - B 1(1.2) 

MATR32 (3,2) - - D1 (2,6) 

IF (1ZZ.EQ.2) GO TO 122 

DO 120 1-1,3 
DO 120 K=l,3 
120 SAVE33 ( I ,K) = MATR33 0.K) 

SAVE3 (1) * NMNY1 
SAVE3 (2) = NMNXY1 
SAVE 3 (3) = NMMY1 
IRR = 0 

CALL LEQT2F (SAVE33. 1 . 3. 3. SAVE 3, 0, WKARE A, I RR) 
M=2 


122 
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N=3 
I A=3 
I RR=0 
I D0=0 

CALL LEQT2F (MATR33,M,N , I A, NATR32 , I DO .WKAREA , I RR) 
DO 125 1-1,3 

IF (I2Z-EQ.2) SAVE30) - ZR 
FNM(I) = SAVE3 (0 
DO 125 L=1 , 2 

125 CD (I , L) ~ MATR32(I,L) 


SC « DSQRT ( ( A1 (5,5) - Al(4,5) * A1 (4,5) / Al(4,4) ) 

C / ( D 1 (6,6) + B1 (2,6) *CD (1 ,2) + B1 (6,6) *CD (2,2) 

C + D1 (2,6) * CD (3.2) ) ) 

GO TO 135 

******************************************************************* 


IN CASE THE LAYERS ARE UNCOUPLED 


********************************************************************* 


130 


c 


DY - -1/ ( A1 (2,2)*A1 (6,6) - A1 (2,6) * A1 (2,6) ) 

CD (1 , 1) - ( A1 (6,6)*A1 (1,2) - A1 (2,6)*A1 (1,6) ) / OY 

CD (2, 1) - ( A1 (2,2)*A1 (1,6) - A1 (2,6) * Al(l,2) ) / OY 

CD (3 , 2) - ( A1 (2,2) *A1 (6,6) + A1 (2,6) *A1 (2,6) ) * 

D 1 (2,6) / D 1 (2,2)/ DY 
CD (1,2) - ZR 

CD (2,2) - ZR 

CD (3, 1) = ZR 


DR - -1 / ( AO (2,2) *A0 (6,6) - AO (2 , 6) *A0 (2 . 6) ) 

WD (1,1) - ( AO(6,6)*AO(l,2) - AO (2 , 6) *A0 (1.6) ) / OR 
WD(2,1) - ( AO (2, 2) *A0 (1,6) - AO (2,6) *A0 (1 ,2) ) / DR 
WD (1 , 2) - ZR 
WD(1,3) “ ZR 
WD (2 ,”2) - ZR 
WD (2 , 3) - ZR 

SSY - DSQRT ( AO ( 5 , 5 ) / DO (6,6) ) 

SSX - DSQRT ( AO (4,4) / DO (2,2) ) 


FNM(1) - (A1 (2,6) *NMNXY1 - A1 (6,6) *NMNY1) / OY 
FNH (2) - (A 1 (2,6) *NMNY1 - A1 (2,2) *NMNXY1 ) / DY 
FNM(3)«(A1 (2,6)*A1 (2,6) - A1 (2,2) *A1 (6,6) ) /DY * NMY1/D1 (2,2) 

SC « DSQRT ( (A1 (5,5)* A1 (6,6)) / (A 1 (6 , 6) *D 1 (6, 6) 

C - (B1 (6,6) *B1 (6.6) ) ) ) 


************************************************************** 

****************************************************************** 


135 Cl « B1(l, 6) + CD (1 , 1) *B1 (2,6) + CD (2, 1) *B1 (6,6) 

C + CD (3. 1) * D 1 (2,6) 

C2 = D1 (6,6) + CO (1,2) *B1 (2.6) + CD (2,2) *B1 (6,6) + 
C CD (3,2) * D 1 (2,6) 


J22 = 00 (2,2) + BO (2,2) * WD(1,2) + BO (2,6) * WO (2,2) 

J66 = DO (6,6) + BO (2 , 6) * WD(1,3) + BO(6,6) * W0(2,3) 

J 26 = DO (2 , 6) + BO (2 , 2) * WD(1,3) + BO(2,6) * WD(2,3) 

BNEG = J22 * AO (5,5) + J66 * AO (4,4) - 2. * J26 * AO (4,5) 

A = J22 * J66 - J26 * J26 



C - AO (4,4) * AO (5,5) - AO (4,5) * AO (14,5) 


SQ •= DSQRT ( (BNEG * BNEG) - I 4 .O * C * A ) 

OIFF1 - DABS ( BNEG + SQ ) 

DIFF2 *= DABS ( BNEG - SQ ) 

IF (DIFF1 .GT.DIFF2) GO TO 140 
SSI = DSQRT ((BNEG + SQ) / 2. / A ) 

GO TO 145 

140 SSI -= DSQRT ( (BNEG-SQ) / 2. / A ) 

145 SS2 « DSQRT ( (BNEG/A) - SSI * SSI ) 

SSY «= DSQRT ( AO (4,4) * J22 ) 

SSX ■= DSQRT ( AO (5,5) * J66 ) 

******************************************************************** 


IF (LIL.EQ.1) WRITE (6,*) ' SI AND S2 * ' , SS 1 , SS2 

IF (LIL.EQ.1) WRITE (6,*)' SX AND SY = \SSX, SSY 

******************************************************************** 


CVNM - ( K66 * (NMNY1 + NMNYO) - K26 * (NMNXY1 + NMNXYO) ) / D 
CUNM - ( K22 * (NMNXY1 + NMNXYO) - K26 * (NMNY1 + NMNYO) ) / D 

CV - STRAIN (LST) / D * ( K26 * Kl6 - K66 * K12 ) + CVNM 
CU «= STRAIN (LST) / D * ( K26 * K12 - K22 * Kl6 ) + CUNM 

****************************************************************** 


NOW FIND SOME OF THE NEEDED CONSTANTS 

FIRST DO LOOP IS TO VARY THE VALUES OF S 


DO 150 I = 1,4 

FORM THE A MATRIX (MATR32) 

MATR33 (1,1) = - ( 
MATR33 (1,2) = - ( 
MATR33(1,3) - ~ ( 
MATR33(2,1) ■= - ( 

MATR33 (2, 2) = - ( 

MATR33 (2,3) - “ ( 

MATR33 (3. 1) - ~ ( 
MATR33 (3,2) - ‘ - ( 
MATR33 (3,3) - - ( 

MATR3 (1) - (F (1,1) 
MATR3 (2) - F (3, 1) 
MATR3(3) - F (2,1) 


AND ITS B (MATR3) 


F ( 3 , 1 ) 

A 

SO) 

* 

SO) 

- 

A 1 (4,5) 

) 

F (2,1) 

* 

SO) 

A 

SO) 

) 



F (4,1) 

* 

s CD 

A 

SO) 

) 



F (3,3) 

* 

SO) 

A 

s(i) 

- 

A l (5,5) 

) 

F (3,2) 

* 

S(l) 

A 

SO) 

) 



F (4,3) 

A 

S(I) 

A 

so) 

) 



F (3.2) 

A 

S(I) 

A 

so) 

) 



F (2,2) 

A 

s(i) 

A 

SO) 

- 

AO (4,4) 

) 

F (4,2) 

* 

SO) 

A 

SO) 

- 

AO (4,5) 

) 

* S (1) 

A 

s(i) 

) 

- A 1 (4,4) 


S(l) 

* S(l) 

- 

A 1 (4,5) 




S ( I ) * S ( I ) 


CALL UP ROUTINE TO FIND THE VALUES OF ALPHA, PHI AND GAMMA 


M=1 

N=3 
I RR=0 
I 00=0 

1A-3 

CALL LEQT2F (MATR33«M,N, I A.MATR3, I DO . WKARE A , I RR) 

ALPHA (I) = MATR3 (1) 

PHI (I) = MATR3 (2) 

GAMMA (I) = MATR3 (3) 

SV(I) = VV11 + ALPHA (l)*VV13 + PHI (l)*VV12 + GAMMA ( I ) *VVl4 
150 SU ( I ) = Ull + ALPHA ( 1 ) *U 1 3 + PHI (l)*U12 + GAMMA(l)*Ul4 


00 155 I - l,A 


155 


NX1 (I)«A1 (1,2)*SV(I) + A1 (1,6)*SU(I) + 
NY1 (l)-Al (2,2)*SV(I) + A1 (2.6) *SU (I) + 
NXY1 (I)*=A1 (2,6) *SV(1) + A1 (6,6)*SU(I) + 
MY1 (I)=B1 (2,2) *SV(I) + B1 (2,6)*SU(I) + 
MXY1 (I)«B1 (2,6) *SV(I) + B1 (6,6) *SU(I) + 
NXO(I) ®A0 (1 ,2) *SV (l) + A0(1,6)*SU(I) + 
GG (1 , I ) = NY 1 (I) 

GG (2, 1) « NXY 1 (I) 

GG (3, 1 ) - MY 1 (I) 

FTH * C2 * SC 

GG (A, I) *= MXY1 (I) + FTH * ALPHA ( 


B1 (1,2) + B1 (1 , 6) *ALPHA ( I ) 
B1 (2,2) + B1 (2,6) *AIPHA (I) 
B1 (2,6) + B1 (6,6) *ALPHA (I) 
D1 (2,2) + D 1 (2,6) *ALPHA (I) 
01 (2,6) + 01 (6,6) *ALPHA (I) 
E 18*PH I (I) + E 1 9*GAMMA ( I ) 


)/ s CD 


A 1 NM - B 1 (2,6)* FNM(l) + B1 (6,6) *FNM (2) + 01(2,6) * FNM (3) 


BG (1) - A1 (1,2) *STRAIN(LST) + CV * A1 (2,2) + A1 (2,6) * CU-NMNY1 

BG (2) « A 1 (1,6) *STRA I N (LST) + CV * A1 (2,6) + A1 (6,6) * CU-NMNXY1 

BG (3) - B1 (1 ,2) ^STRAIN (LST) + B1 (2,2) * CV + B1 (2,6) * CU-NMMY1 

BG (A) - B1 (1,6) *STRA I N (LST) + B1 (2,6) * CV + 61 (6,6) * CU-A1NM 

C - Cl * STRAIN (LST) 

BG 1= BG (1) 

BG2 - BG (2) 


M=1 

n=a 

I A=A 
I DD=0 
1 RR=0 

CALL LEQT2F (GG ,M, N, I A, BG , I 00 .WKAREA, 1 RR) 
*************>v***********5ftr***yc*********^*****jv*Ajv***/f**************yt 

TVNM = S1NM - FNM(l) + HI / 2.0 * FNM (3) 

TUNM = S2NM - FNM (2) 

THETV * -CD (1,1) + WD (1 , 1) + H 1/2 .0*CD (3, 1) 

THETU - -CO (2,1) + WD (2 , 1 ) 

THETV = THETV + TVNM 
THETU ■= THETU + TUNM 

IF (LIL.EQ.1) WRITE (6,215) THETV, THETU 
IF (LIL.EQ.1) WRITE (6,216) SMNY, SMNXY 

THE STEPS USED TO FIND THE TOTAL ENERGY RELEASE FROM USING A PURE 
EXTENSION ANALYSIS FOR THE HYGRALTHERMAL EFFECTS. (SIMILAR TO 
WHITNEY'S). ANALYSIS 1$ CARRIED OUT ON A PLY BY PLY BASIS 


ZV - ( K26 * K 1 6 - K66 * K12 ) / 0 
ZU = ( K26 * K12 - K22 * Kl6 ) /D 
C 

DO 162 LL ■= l.TPLY 
C 

c 

EX (LL) = THICK(LL) * ( Q(1,1,LL) + ZV*Q(1,2,LL) + 

C ' ZU*Q(1,6,LL) ) 

Tl(LL) = NXNM(LL) - CVNM*TH I CK (LL) *Q (1 , 2, LL) - 
CUNM * THICK (LL) * Q(1 ,6,LL) 

EXX (LL) = Q(1,1,LL)*THICK(LL) + Q (1 , 2 , LL) *TH I CK (LL) *CD (1 , 1) 


C 



c 


+ Q(1,6,LL)*THICK(LL)*CD(2,1) + B12(LL)*CO(3,l) 

T1 2 (LL) ■= NXNM(LL) - FNM (1) *Q (1 , 2, LL) *TH I CK (LL) - 

C FNM (2) *Q (1 ,6, LL) *TH I CK (LL) - FNM (3) *B 1 2 (LL) 

EX3(LL) = Q(1,1,LL)*THICK(LL) + WO (1 , 1) *Q (1 ,2, LL) *THI CK (LL) 
C + WD(2,1)*Q(1,6,LL)*THICK(LL) 

T1 3 (LL) = NXNM(LL) - Q (1 , 2 , LL) *TH I CK (LL) *S1NM - Q(1,6,LL) 

C * TH I CK (LL) * S2NM 

IF (I ZZ.EQ.2) NMSTO(LL) « 0.0 
IF (I ZZ.EQ.2) NMST2 (LL) = 0.0 

IF (I ZZ.EQ.2) NMST3 (LL) = 0.0 

IF (I ZZ.EQ.2) GO TO 162 


NMSTO(LL) = T1 (LL) / EX (LL) 

NMST2 (LL) = T1 2 (LL) / EXX (LL) 

NMST3 (LL) = T1 3 (LL) / EX3(LL) 

162 CONTINUE 

WRITE (6,*)' JMM, NMSTO, 2. 3 OF ALL PLYS ',( NMSTO (JP) , 
C NMST2 (JP) ,NMST3(JP) , ' — ' ,JP«1,TPLY) 

WRITE (6,*)’ EX, EXX EX3 OF ALL PLYS *,( EX (JP) , 

C EXX(JP) ,EX3(JP) 1 , JP=1 ,TPLY) 


163 


TNC - 0.0 
EXNC = 0.0 
TSTAR =0.0 
ESTAR = 0.0 
DO 163 LK = l.TPLY 
TNC = TNC + T1 (LK) 

EXNC = EXNC + EX (LK) 

IF (LK.LE.NPLYO) TSTAR = TSTAR + T13(LK) 
IF (LK.GT.NPLYO) TSTAR = TSTAR + T12(LK) 


IF (LK.LE.NPLYO) ESTAR 
IF (LK.GT.NPLYO) ESTAR 


ESTAR + EX3(LK) 
ESTAR + EXX (LK) 


IF (l ZZ.EQ.2) TNMST = 0.0 
IF (I ZZ.EQ.2) GO TO 89 


TNMST = ( TNC - (TNC-TSTAR) *2*AL/WI DTH ) / ( EXNC - 
C (EXNC-ESTAR) *2*AL/WI DTH ) 

89 WRITE (6,*)' TNMST EQUALS ' .TNMST 


DO 164 LL = l.TPLY 

IF (LL.LE.NPLYO) WWC = (EX3(LL) * 

C (STRAIN(LST) + TNMST) - T13(LL) ) 

C * ( STRAIN(LST) - NMST3 (LL) + TNMST ) 

IF (LL.GT.NPLYO) WWC = (EXX(LL) * 

C (STRAIN(LST) + TNMST) -T12(LL) ) 

C * ( STRAIN(LST) - NMST2 (LL) + TNMST ) 

WWO = ( EX (LL) * (STRAIN(LST) + TNMST) - T1 (LL) ) 
C * ( STRAIN(LST) - NMSTO (LL) + TNMST ) 


164 GLC(JMM) = GLC(JMM) + WWO - WWC 


GLC(JMM) = GLC(JMM) / 2.0 

************************** * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 



THIS IS TO CALCULATE THE INTERLAMINAR SHEAR STRESSES. 

************************************************************* ******** 
UNCL - (WIDTH / 2.0)- AL 
DO 180 JX - 80,100 
JY «( 1.0 - JX /100.0) * UNCL 

SIGX(JMM.JX) - 0 
SI GY (JMM, JX) « 0 
DO 180 JS = 1,4 

SIGX(JMM.JX) «' BG(JS) * S (JS) * DEXP ( -S (JS) * JY ) 

C * NXY 1 (JS) + SIGX (JMM, JX) 

S I GY (JMM, JX) = BG (JS) * S (JS) * DEXP ( -S (JS) * JY ) 

C * NYl(JS) + SI GY (JMM, JX) 

180 CONTINUE 


THE PROGRAM CONTINUES AND FINDS THE VARIOUS STRAIN ENERGY RELEASE 
COMPONENTS 


IF (COUPL.EQ.l) GO TO I 65 

THIS IS FOR A SYSTEM THAT IS COUPLED, THE CRACK LENGTH IS 

DEL - S(4) * S (2) * ATH * ATH 
DEL - DEL * DEL * 0.6144 
GO TO 170 

THIS IS FOR AN UNCOUPLED SYSTEM 

I 65 SSW = .65 * ( S(l) + S( 2 ) + S ( 3 ) ) 

DEL - 18.7 * S(4) * SSW * ATH * ATH 
DEL = DEL * DEL / 571-00 
170 DEL - 135-7 * DEL * ATH 

C IF (LIL.EQ.1) WRITE (6,21 1) DEL 

FY - ZR 
FX - ZR 


DO 175 JP=1 ,4 

CONY- NY1 (JP) *BG (JP) * ( DEXP ( -S (JP) * DEL ) - 1 ) 
CONXY- NXY 1 (JP) * BG (JP) * ( DEXP ( -S (JP) *DEL)- 1 ) 
CCONXY - CCONXY + CONXY/S (JP) 

175 CCONY - CCONY + CONY/S (JP) 

*************************************************************** 

*************************************************************** 

FY - BG1 + CCONY / DEL 
FX - BG2 + CCONXY / DEL 

GII(IZZ.JM) = FY / 2.0 * THETV *STRAIN(LST) 

Gill (IZZ.JM) = FX / 2.0 * THETU * STRAIN (LST) 

DIFFG- Gl I (IZZ, JM) - G I I I (I ZZ, JM) 

CON = 2 

IF (DIFFG. GT. REAL (GLC (JMM))) CON-1 

IF (DIFFG. GT. REAL (GLC (JMM))) DEL = DEL * -9 

IF (DIFFG. GT. REAL (GLC (JMM))) WRITE( 6 ,*)' IT EXPLODES 


Gl (IZZ.JM) = GLC (JMM) -Gl I (IZZ, JM) -Gill (IZZ.JM) 



200 


CONTINUE 


RESULTS ARE PRINTED FOR EACH RUN. 


300 CONTINUE 

WRITE (6,266) STRAIN(LST) 

WRITE (6,267) 

WRITE ( 6 , 269 ) GLC(O) , G I (2.1) , G I l (2.1). GMI(2.1), 

C G I (2,1) /GLC (0) 

DO 350 l-l.MMC +1 

WRITE (6,268) CMOIST(I) ,GLC (I) , G I (1 , I ) , G I I (1 , I ) , 

C Gl I I (1, I) , Gl (1 , I) /GLC (I) 

350 CONTINUE 

WRITE (6,287) 

DO 36 O NS - 80,90,2 

360 WRITE (6,285)NS/100., S I GX (0 , NS) , ( S I GX (KL , NS) , KL«1 , 22 , 4) 

DO 365 NS « 91,100 

365 WR I TE ( 6 , 285) NS/ 100 . , S I GX ( 0 , NS) , ( S I GX (KL , NS) , KL-1 , 22 , 4) 

WRITE (6,286) 

DO 370 NS - 80,90,2 

370 WRITE(6,285) NS/100. ,SIGY(0, NS) ,( SIGY(KL, NS) , KL-1, 22, 4) 

DO 375 NS « 91.100 

375 WRITE(6,285) NS/ 100 . ,SIGY(0, NS) ,( SIGY(KL, NS) , KL«1,22,4) 


400 CONTINUE 

201 FORMAT (//,' THE WIOTH OF THE LAMINATE IS \F 8 . 5 ) 

287 FORMAT (/////, ' THESE ARE THE IN-PLANE INTERLAMINAR SHEAR 
C 'STRESSES -- SIGMA XY ',/,' THEY ARE FOUND AT VARIOUS.', 

C '• MOISTURE CONTENTS •',//,' Y LOCAT I ON • , 7X, ' MECH 0NLY',8X, 

C 'H-0.0' , 10X, 'H-0.2' ,10X, 'H-0.4' , 10X, 

C 'H-0.6' , 10X, 'H-0.8' , 10X, 'H-l .0' ,//) 

285 FORMAT (3X,F7.2,4X,7F15*8) 

202 FORMATC THE NUMBER OF LAMINATES ABOVE AND BELOW THE CRACK IS' 

C . 13 . 5 x. 13 ) 

204 FORMAT (///,' THE PLYS ARE INPUTTED FROM BOTTOM TO TOP',/, 

. C 'BUT THE PLY CHARACTERISTICS FROM TOP TO BOTTOM ARE ') 

206 FORMAT (//, ' FOR PLY ',15,' THE SUBLAMINATE HAS THESE PROPERTIES') 

205 FORMAT (//,' WITH THIS LAYUP, THE PLYS ARE COUPLED ',//) 

210 FORMAT (//,' WITH THIS LAYUP, THE PLYS ARE DECOUPLED ',//) 

286 F ORMAT (//, ' THESE ARE THE OUT-OF-PLANE INTERLAMINAR SHEAR ', 

C 'STRESSES — SIGMA YZ ',/,' THEY ARE FOUND AT VARIOUS', 

C ' MOISTURE CONTENTS ',//,' Y LOCAT I ON ' ,8X, ' MECH 0NLY',8X, 

C * H=0 . 0 ' , 1 OX , ' H=0 . 2 ' , 1 OX , ' H=0 . 4 ' , 10X, 

C 'H-0.6' , 10X, 'H-0.8' ,10X, 'H-l .0' ,//) 

208 FORMATC El AND E2 ARE (MSI) ' ,F8.4, 10X.F8.4) 

289 FORMAT (//,' THE LAMINA PLY CHARACTERISTICS INITIALLY ARE ',/) 

288 FORMAT (/,' S1NM AND S2NM ARE EQUAL TO ' , F 1 4 . 10, 4X , F 1 4 . 10) 

211 FORMAT (//,' THE CRACK LENGTH STEP SIZE IS ' , F 1 2 . 8) 

266 FORMAT ('O' T , ' THE STRAIN IS EQUAL TO ' , F 1 2 . 7 , /, 

C ' THE VALUES OF GT, Gl, Gil, AND GUI ARE IN IN-LB/IN/IN 
C ') 

267 FORMAT (/,3X, '% CMOIST' ,8X, 'GGG (WHITNEY) ' ,6X, 'Gl ' ,9X, 

C ' Gl I ' ,8X, 'Gl I I ' ,6X, ' Gl/G (W-T) ' ,//) 

269 FORMAT (/, ' MECH. ONLY ' , 3X,F 12.9.4 ( 2 X,F 1 1 . 7 ) ,/) 

268 FORMAT(5X,F8.3,3X,F12.9.4(2X,F11.7) ) 

215 FORMAT (////, ' THETA V. IS ' , F 15 . 10 , ' THETA U IS ' , F 24 . 1 9) 

216 FORMAT (//, ' NY IS '.F23.ll,' NXY IS ' , F23 - 1 8) 



207 FORMAT (/,' THE THICKNESS AND THETA VALUES ARE ' , F9- 6 .5X. F8 . 3) 
209 FORMAT ( 1 THE POISSON RATIO (1.2) IS ' .F10.5) 

211* FORMAT ( ' G OF (1-2), AND (3~1) ARE -MSI \ 2(F9-A,2X)) 

217 F0RMAT('0' ,////, 8X, 'THE FOUR CHARACTERISTIC VALUES ASSOCIATED' 

C ,/,8X, ' WITH THE 8 DEGREE POLYNOMIAL FOR THE COUPLED CASE ARE') 

218 FORMAT (//,5X. ' THETA ' ,£>X, ' NY ' ,8X, ' NXY ' ,//) 

219 FORMAT (///,' THE S VALUES OF THE MEMBRANE ARE ' , F 1 5 -5. 3X, F 15 -5) 

220 FORMAT(/,LX,F9.l4,3X,F9.2,3X,F9.2) 

221 FORMAT (/, ' S OF ',12,' IS EQUAL TO \F20.10) 

223 FORMAT (/, ' THE UNCOUPLED SY (1,2) VALUES ARE ' , F 15 -5 . 3X, F 1 5 -5) 

22A FORMAT (/,' THE UNCOUPLED SX (1,2) VALUES ARE * ,F 15-5.3X.F 15-5) 

231 FORMAT (//,' THE STRAIN IS EQUAL TO '.F12.8,/, 

C ' THE CHANGE IN TEMPERATURE IS '.F12.5,/, 

C ' THE COEFFICIENTS SWELLING DUE TO MOISTURE ARE ' , 2 (2X, F 1 2 .8) 

C THE COEFFICIENTS OF THERMAL EXPANSION ARE ' ,2 (2X,F I 5 . 9 ) ) 


232 FORMAT (/,' THE MOISTURE COEFF I CENT IS ' . F 1 5 . 8 ) 

233 FORMAT (/,' THE MOISTURE COEFFICIENT VARIES FROM 0 TO 1.2 ') 



THE LAMINA PLY CHARACTERISTICS INITIALLY ARE 
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THE WIDTH OF THE LAMINATE IS 1.51200 
THE NUMBER OF LAMINATES ABOVE AND BELOW THE CRACK IS 3 1 


FOR PLY 1 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .0051*00 35.000 

El AND E2 ARE (MSI) 18.7000 1.2300 

THE POISSON RATIO (1.2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


FOR PLY 2 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE . 005^00 -35-000 

El AND E2 ARE (MSI) 18.7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


FOR PLY 3 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .0051*00 
El AND E 2 ARE (MSI) I 8 . 7 OOO 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1)' ARE -MSI .8320 .8320 


.000 

1.2300 


FOR PLY A THE SUBLAMINATE HAS THESE PROPERT 
THE THICKNESS AND THETA VALUES ARE .0051*00 
El AND E2 ARE (MSI) 18.7000 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 


IES 

90.000 

1.2300 


.8320 


THE STRAIN IS EQUAL TO .00251*000 
THE CHANGE IN TEMPERATURE IS -280.00000 

THE COEFFICIENTS SWELLING DUE TO MOISTURE ARE .00000000 

THE COEFFICIENTS OF THERMAL EXPANSION ARE -.000000230 

THE MOISTURE COEFFICIENT VARIES FROM 0 TO 1.2 


WITH THIS LAYUP, THE PLYS ARE COUPLED 


S 

S 

S 

S 


THE FOUR CHARACTERISTIC- VALUES ASSOCIATED 
WITH THE 8 DEGREE POLYNOMIAL FOR THE COUPLED CASE ARE 


OF 

1 

IS 

EQUAL TO 

1 * 07 . 057368271 * 1 * 

OF 

2 

IS 

EQUAL TO 

11 * 1.1 1977801*18 

OF 

3 

IS 

EQUAL TO 

1 16.7332723860 

OF 

1 * 

IS 

EQUAL TO 

55.551*1*207729 


.00556000 

. 000011*900 


THE UNCOUPLED SX (1,2) VALUES ARE 


392.221*78 


106.21737 


THE UNCOUPLED SY (1,2) VALUES ARE 


13^.68589 


56 . 2575 ** 


THE S VALUES OF THE MEMBRANE ARE 1 77 - 389^5 66.26466 


S1NM AND S2NM ARE EQUAL TO -.0000157292 .0000000000 

SI AND S2 = 1 34.932471163957803943120257 641.500299099584182787943089 
SX AND SY = 7.OO3582O8142O9O671392409983 33.2966554398959572557143052 


THETA V IS .3932559441 THETA U IS -.0981017742335529623 
NY IS -38.32041690358 NXY IS -61.959633961712929972 


THE STRAIN IS 

EQUAL TO .0025400 

THE VALUES OF 

GT, Gl, Gil, 

AND GUI 

% CMOIST 

GGG (WHITNEY) 

Gl 

MECH. ONLY 

.101653935 

.0670191 

.000 

•522939955 

.4084740 

.050 

.488576635 

•3794622 

.100 

.455156902 

.3513864 

.150 

.422680755 

.3242463 

.200 

. 39"S 148194 

.2980421 

.250 

.360559220 

.2727737 

.300 

.330913832 

.2484412 

• 350 

.302212031 

.2250445 

.400 

.274453816 

.2025836 

.450 

.247639187 

.1810586 

.500 

.221768145 

.1604694 

.550 

.196840689 

.1408161 

.600 

.172856820 

.1220986 

.650 

.149816537 

.1043169 

• 700 

.127719840 

.0874710 

.750 

. 106566730 

.0715610 

.800 

.086357206 

.0565869 

.850 

.067091269 

.0425486 

.900 

.048768918 

.0294461 

• 950 

.031390154 

.0172794 

1 .000 

.014954976 

.0060486 

1.050 

-.000536616 

-.0042464 

1.100 

-.015084621 

-.0136055 

1.150 

-.028689040 

-.0220288 

1 .200 

-.041349872 

-.0295163 


ARE IN IN-LB/IN/IN 


Gl 1 

Gl 1 1 

Gl/G (W-T) 

.0346174 

.0000174 

.6592868 

.1144151 

.0000509 

.7811106 

. 1090658 

.0000486 

•7766688 

.1037241 

.0000464 

.7720115 

.0983903 

.0000442 

.7671187 

.0930641 

.0000420 

.7619672 

.0877458 

.0000397 

.7565296 

.0824351 

.0000375 

.7507731 

.0771323 

.0000353 

.7446576 

.0718371 

.0000331 

.7381338 

.0665497 

.0000308 

.7311387 

.0612701 

.0000286 

.7235909 

.0559982 

.0000264 

.7153809 

.0507341 

.0000242 

.7063566 

»0454777 

.0000220 

.6962975 

.0402291 

.0000197 

.6848665 

.0349882 

.0000175 

.6715139 

.0297550 

.0000153 

.6552653 

.0245296 

.0000131 

.6341891 

.0193120 

.0000109 

.6037875 

.0141021 

.0000087 

.5504724 

.0088999 

.0000064 

.4044538 

.0037055 

.0000042 

7.9132595 

-.0014811 

.0000020 

.9019465 

-.0066600 

-.0000002 

.7678481 

-.0118312 

-.0000024 

.7138182 



THESE ARE THE IN-PLANE INTERLAMINAR SHEAR STRESSES -- SIGMA XZ 
•l HEY ARE FOUND AT VARIOUS MOISTURE CONTENTS 
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THE LAMINA PLY CHARACTERISTICS INITIALLY ARE 


THE WIDTH OF THE LAMINATE IS 1.51200 
THE NUMBER OF LAMINATES ABOVE AND BELOW THE CRACK IS 3 1 

THE PLYS ARE INPUTTED FROM BOTTOM TO TOP 
BUT THE PLY CHARACTERISTICS FROM TOP TO BOTTOM ARE 


FOR PLY 1 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .0051*00 35.000 

El AND E2 ARE (MSI) 18-7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2). AND (3-1) ARE -MSI .8320 .8320 


FOR PLY 2 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE . 005^00 

El AND E2 ARE (MSI) 18-7000 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


.000 

1.2300 


FOR PLY 3 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .005^00 -35. 

El AND E2 ARE (MSI) 18.7000 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


000 

1.2300 


FOR PLY It THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .005^00 90. 

El AND E2 ARE (MSI) • 18.7000 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


000 

1.2300 


THE STRAIN IS EQUAL TO .00251*000 
THE CHANGE IN TEMPERATURE IS -280.00000 

THE COEFFICIENTS SWELLING DUE TO MOISTURE ARE .00000000 

THE COEFFICIENTS OF THERMAL EXPANSION ARE -.000000230 

THE MOISTURE COEFFICIENT VARIES FROM 0 TO 1.2 


WITH THIS LAYUP, THE PLYS ARE COUPLED 


THE FOUR CHARACTERISTIC VALUES ASSOCIATED 
WITH THE 8 DEGREE POLYNOMIAL FOR THE COUPLED CASE ARE 


OF 

1 

IS 

EQUAL 

TO 

360. 71 621*235^3 

OF 

2 

IS 

EQUAL 

TO 

136.396160LL92 

OF 

3 

IS 

EQUAL 

TO 

113-9856581*772 

OF 

*4 

IS 

EQUAL 

TO' 

5b. 080 1738692 


.00556000 

. 0000 U 900 



THE UNCOUPLED 

SX (1,2) VALUES 

ARE 

352.6287*4 

86.89277 


THE UNCOUPLED 

SY (1,2) VALUES 

ARE 

126. *45615 

59.53405 

63 

THE S VALUES OF 

THE MEMBRANE ARE 193-07807 

70.65819 


S1NM AND S2NM 

ARE EQUAL TO 

■.0000157292 .0000000000 


SI AND S2 

= 1314.932*4711639578039^3120257 614 1 .50029909958*4 1827879*43089 

SX AND SY 

■= 7.003582081142090671392*409983 33.296655*43989595725571*43052 

THETA V IS 

• 978*402*4562 THETA U IS 

.00000000000000097 1 8 


NY IS 

-38.320*4 1690358 

NXY IS - 

■61.959633961712929972 


THE STRAIN IS 

EQUAL TO .0025*400 




THE VALUES OF 

GT, Gl, Gil, 

AND Gill 

ARE IN IN- 

LB/IN/IN 


% CMO 1 ST 

GGG (WHITNEY) 

Gl 

Gl 1 

Gill 

Gl/G (W-T) 

MECH. ONLY 

.09*4207177 

.0078275 

.0863797 

.0000000 

.0830882 

.000 

.510289070 

.2262277 

.284061*4 

.0000000 

.4433324 

.050 

.*476*405835 

.2055316 

.2708742 

.0000000 

.4314212 

. TOO 

.*♦143*4*46720 

. 1857*497 

.2576970 

.0000000 

.4188771 

.150 

.*411*411727 

. 1668819 

.2445298 

.0000000 

.4056324 

.200 

.380300856 

. 1*48928*4 

.2313725 

.0000000 

.3916067 

.250 

.3501 1*4106 

.1318890 

.2182251 

.0000000 

.3767028 

.300 

.320851*477 

•1157637 

.2050878 

.0000000 

.3608016 

• 350 

.292512969 

.1005526 

.1919603 

.0000000 

.3437545 

. 400 

.265098582 

.0862557 

.1788428 

.0000000 

.3253723 

.1*50 

.238608317 

.0728730 

•1657353 

.0000000 

.3054085 

.500 

.2130*42173 

.060*40*4*4 

•1526377 

.0000000 

.2835328 

.550 

. 188*400151 

.01488500 

.1395501 

.0000000 

.2592888 

.600 

. 16*4682250 

.0382098 

.1264724 

.0000000 

.2320215 

.650 

. 1*41888*470 

.028*4838 

.1134047 

.0000000 

.2007475 

.700 

.12001881 l 

.0196719 

.1003469 

.0000000 

.1639065 

.750 

.09907327*4 

.01177*41 

.0872991 

.0000000 

.1188427 

.800 

.079051858 

.00*47906 

.0742613 

.0000000 

.0606005 

.850 

.05995*4563 

-.0012788 

.0612334 

.0000000 

-.0213295 

.900 

.0*4 1 781389 

-.006*43*40 

.0482154 

.0000000 

-.1539925 

.950 

.02*4532337 

-.0106751 

.0352074 

.0000000 

-.435142° 

1 .000 

.008207*406 

-.01*40020 

.0222094 

.0000000 

-1.7060144 

1.050 

-.007193*403 

-.016*41*47 

.0092213 

.0000000 

2.2819054 

1 . 100 

-.021670092 

-.0179132 

-.0037569 

.0000000 

.8266329 

1.150 

-.035222659 

-.018*4976 

-.0167251 

.0000000 

.5251615 

1 .200 

-.0*47851 10*4 

-.0181678 

-.0296833 

.0000000 

.3796734 


THESE ARE THE IN-PLANE INTERLAMINAR SHEAR STRESSES -- SIGMA XZ 
THEY ARE FOUND AT VARIOUS MOISTURE CONTENTS 
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THE LAMINA PLY CHARACTERISTICS INITIALLY ARE 


85 


THE WIDTH OF THE LAMINATE IS 1.51200 
THE NUMBER OF LAMINATES ABOVE AND BELOW THE CRACK IS 2 


FOR PLY 1 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .005400 30.000 

El AND E2 ARE (MSI) 18.7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


FOR PLY 2 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .0051*00 -60.000 

El AND E2 ARE (MSI) 18.7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


FOR PLY 3 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .005400 75-000 

El AND E2 ARE (MSI) 18.7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


FOR PLY 4 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .005400 -15-000 

El AND E2 ARE (MSI) 18.7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


THE STRAIN IS EQUAL TO .00251*000 
THE CHANGE IN TEMPERATURE IS -280.00000 

THE COEFFICIENTS SWELLING DUE TO MOISTURE ARE .00000000 .00556000 

THE COEFFICIENTS OF THERMAL EXPANSION ARE -.000000230 .000014900 

THE MOISTURE COEFFICIENT VARIES FROM 0 TO 1.2 


WITH THIS LAYUP, THE PLYS ARE COUPLED 


THE FOUR CHARACTERISTIC VALUES ASSOCIATED 
WITH THE 8 DEGREE -POLYNOMIAL FOR THE COUPLED CASE ARE 

S OF IIS EQUAL TO 202 . 3666962066 

S OF 2 IS EQUAL TO 150.2447231*209 

S OF 3 IS EQUAL TO 96.3990023366 

S OF 4 IS EQUAL TO 72.1855545293 


THE UNCOUPLED S>: (1.2) VALUES ARE 


178.08581 


9C.73315 



THE UNCOUPLED SY (1,2) VALUES ARE 


137.79321. 


77. 57481 


THE S VALUES OF THE MEMBRANE ARE 107-30161. 75-55664 

S1NM AND S2NM ARE EQUAL TO -.0013126127 .0012407192 

SI AND S2 « li«5 -7 12738391*363 11*5201 23A 301* 261.71259626539561223751*2175 
SX AND SY ■= 35.5756851 1671*391 1519011*7879 60.958551*9677311*312319553312 


THETA V IS .2402517909 THETA U IS 2 . 07 31*381+470 2 1 801786 1 
NY IS -35.1*8501616832 NXY IS -6l.l*6l8509109l*0726086 


THE STRAIN IS EQUAL TO .00251*00 


THE VALUES OF 

GT, Gl, Gil, 

AND Gill 

ARE IN IN- 

LB/IN/IN 


S CMOIST 

GGG (WHITNEY) 

Gl 

G 1 1 

Gl 1 1 

G 1 /G (W-T) 

MECH. ONLY 

.171*065036 

.0136731* 

.0093651 

.1510265 

.0785535 

.000 

.28859921*1* 

.0062276 

.0365562 

• 21*58155 

.0215787 

.050 

.277169512 

.0030909 

.031*71*33 

.2393353 

.011 1516 

.100 

.266289511* 

.0005010 

.0329323 

.2328563 

.0018813 

.150 

.255959252 

-.00151*22 

.0311230 

.2263784 

-.0060250 

.200 

.21.6178721* 

-.0030385 

.0293158 

.2199016 

-.0123427 

.250 

.23691*7932 

-.0039880 

.0275100 

.2131*260 

- .0168309 

.300 

.228266871* 

-.001*3908 

.0257062 

.2069515 

-.0192354 

• 350 

•220135552 

-.001*21*68 

.023901*2 

.2004781 

-.0192916 

.1*00 

.212553961* 

-.0035559 

.022101*0 

.1940059 ‘ 

-.0167296 

.1*50 

.205522112 

-.0023183 

.0203057 

.1875347 

-.0112801 

.500 

.199039991* 

-.0005339 

.0185092 

.1810647 

-.0026824 

• 550 

.193107612 

.0017973 

.016711*1* 

.1745959 

.0093072 

.600 

. 187721*965 

.001*6753 

.011*9215 

.1681281 

.0249050 

.650 

. 182892052 

.0081001 

.0131301* 

.1616615 

.0442888 

.700 

. 178608875 

.0120716 

.01131*12 

.1551961 

.0675870 

• 750 

.17i*875i*32 

.0165900 

.0095537 

.1487317 

.0948676 

.800 

.171691725 

.0216552 

.0077681 

. 1422685 

.1261282 

.850 

.169057753 

.0272671 

.005981*2 

.1358064 

. 1612888 

• 900 

.166973515 

.0331*259 

.001*2022 

.12931*51* . 

.2001866 

• 950 

.1651*39013 

.01*01311* 

.0021*220 

.1228856 

.2425751 

1 .000 

.161*1*51*21*5 

.01*73837 

.0006L36 

.1164269 

.2881270 

1 .050 

.161*019213 

.0551828 

-.0011329 

.1099693 

.3364412 

1 . 100 

. 161*133916 

.0635287 

-.0029077 

.1035129 

.3870542 

1.150 

.161*798353 

.0721*211* 

-.001*6806 

.0970576 

.4394548 

1.200 

.166012526 

.0818609 

-.0061*517 

.0906034 

.4931008 



THESE ARE THE IN-PLANE INTERLAMINAR SHEAR STRESSES -- SIGMA XZ 
THEY ARE FOUND AT VARIOUS MOISTURE CONTENTS 
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THE LAMINA PLY CHARACTERISTICS INITIALLY ARE 


THE WIDTH OF THE LAMINATE IS 1-51200 88 

THE NUMBER OF LAMINATES ABOVE AND BELOW THE CRACK IS 2 2 


FOR PLY 1 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .0051*00 -35-000 

El AND E2 ARE (MSI) 18-7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


FOR PLY 2 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .005400 55-000 

El AND E2 ARE (MSI) 1 8 . 7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


FOR PLY 3 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .005400 10.000 

El AND E2 ARE (MSI) 1 8. 7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


FOR PLY 4 THE SUBLAMINATE HAS THESE PROPERTIES 
THE THICKNESS AND THETA VALUES ARE .005400 -80.000 

El AND E2 ARE (MSI) 18.7000 1.2300 

THE POISSON RATIO (1,2) IS .29200 

G OF (1-2), AND (3-1) ARE -MSI .8320 .8320 


THE STRAIN IS EQUAL TO .00254000 
THE CHANGE IN TEMPERATURE IS -280. 00000 

THE COEFFICIENTS SWELLING DUE TO MOISTURE ARE 
THE COEFFICIENTS OF THERMAL EXPANSION ARE 


.00000000 

-.000000230 


.00556000 

.000014900 


THE MOISTURE COEFFICIENT VARIES FROM 0 TO 1.2 


WITH THIS LAYUP, THE PLYS ARE COUPLED 


THE FOUR CHARACTERISTIC VALUES ASSOCIATED 
WITH THE 8 DEGREE POLYNOMIAL FOR THE COUPLED CASE ARE 


S OF 

1 

IS EQUAL 

TO 

233.9388572236 


S OF 

2 

IS EQUAL 

TO 

156.9014619788 


S OF 

3 

IS EQUAL 

TO 

115.2992565645 


S OF 

4 

IS EQUAL 

TO 

48.0575100718 


THE 

UNCOUPLED SX 

(1.2) 

VALUES ARE 186.41344 

96.42124 

THE 

uncouple:- S'- 


VALUES ARE 128.78572 

49.37034 



THE S VALUES OF THE MEMBRANE ARE 


1 19.36060 


57.91987 


89 


S1NM ANO S2NM ARE EQUAL TO - .0007221 149 -.0007139234 

SI AND S2 - H3.69l863395728751259l*52089 285.261.127101.852281 15901*5526 
SX AND SY - 32. 1436979546791447820272102 62.201.9978939766992050711339 


THETA V IS .1*31.7536621 THETA U IS -1 .731*44 151.08201 789576 
NY IS -51.. 366 1988991.6 NXY IS 1.5.6186571.1.5050692111 


THE STRAIN IS EQUAL TO .00251.00 


THE VALUES OF 

GT, Gl, Gil, 

AND Gill 

ARE 4 N IN- 

•LB/IN/IN 


% CMOIST 

GGG (WHITNEY) 

Gl 

G 1 1 

Gill 

G 1 /G (W-T) 

MECH. ONLY 

.131525889 

.0207693 

.0219689 

.0887877 

.1579102 

.000 

.391.637005 

.1683076 

.0866210 

.1397081* 

.1*261*871 

.050 

.36931.7077 

.1508205 

.0823111* 

.1362152 

.40831*35 

.100 

.3^5181*223 

.131*1*556 

.0780061 

.1327225 

.3895185 

. 150 

.32211*81*1*3 

.1192131 

.0737050 

. 1292301* 

.3700561* 

.200 

.300239735 

. 1050929 

.0691*082 

.1257387 

.3500298 

.250 

.2791*58102 

.092091*9 

.0651156 

. 12221*75 

.3295483 

.300 

.2598035lt2 

.0802193 

.0608271* 

.1187569 

.3087691 

•350 

.21*1276055 

.0691.660 

.05651*31* 

.1152667 

.2879108 

.1.00 

.22387561*2 

.0598350 

.0522636 

•1117771 

.2672688 

.1*50 

.207602302 

.0513262 

.01*79882 

. 1082879 

.2472335 

.500 

.1921*56036 

.01*39398 

.01*37170 

. 101*7992 

.22831 10 

.550 

. 1781*3681*3 

.0376757 

.0391*500 

.1013111 

.21 1 1431 

.600 

. l655l*l*721* 

.0325339 

.0351871* 

.0978235 

.1965261* 

.650 

.153779678 

.028511*1* 

.0309290 

.091.3363 

.1854237 

.700 

. 1 1*3 1 1+ 1 706 

.0256172 

.026671*8 

.09081*97 

.1789638 

.750 

.133630807 

.02381*23 

.0221*250 

.0873635 

.1784190 

.800 

.12521*6982 

.0231897 

.0181791* 

.0838779 

.1851516 

.850 

.117990230 

.0236591* 

.0139381 

.0803928 

.2005198 

.900 

. 1 11860552 

.0252511* 

.0097010 

.0769081 

.2257398 

•950 

. 10685791*7 

.0279657 

.0051*683 

.073421.0 

.2617089 

1 .000 

. 1029821*16 

.0318023 

.0012397 

.06991*01* 

.3088128 

1 .050 

.100233958 

.0367612 

-.002981*5 

.0661*573 

•3667539 

1.100 

.098612573 

• Ol.281.2l. 

-.007201*5 

.0629747 

.4344517 

1.150 

.0981 18262 

.05001*59 

-.01 14202 

.0594925 

.5100571 

1 .200 

.098751025 

.0583717 

-.0156316 

.0560109 

.5910999 



THESE ARE THE IN-PLANE INTERLAMINAR SHEAR STRESSES -- SIGMA X2 
THEY ARE FOUND AT VARIOUS MOISTURE CONTENTS 
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